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SUMMARY OF METHODS FOR CALCULATING DYNAMIC LATERAL STABILITY AND RESPONSE 
AND FOR ESTIMATING LATERAL STABILITY DERIVATIVES 1 

Bv John P. Campbell and Marion O. McKinney 


SUMMARY 

A summary of methods for making dynamic lateral stability 
and response calculations and for estimating the aerodynamic 
stability derivative# required for use in these calculations is pre- 
sented. The processes of performing calculations of the time 
histories of lateral motions , oj the period and (lam piny of these 
motions , and of the lateral stability boundaries are presented as 
a series of simple straightforward steps . Existing methods for 
estimating the stability deri rati res are summarized and , in 
some cases, simple new em pineal formulas are presented. Ref- 
erence is also made to reports presenting experimental data 
that should be useful in making estimates of the derivatives. 
Detailed estimation methods are presented for low-subsonic- 
speed conditions but only a brief discussion and a list of refer- 
ences are given for transonic- and supersonic-speed conditions. 

INTRODUCTION 

Dynamic lateral stability has not received widespread 
attention in the past because it has not generally been a 
serious problem in the design of airplanes. Consideration 
of dynamic lateral stability has recently become more im- 
portant, however, because current design trends toward the 
use of low aspect ratio, sweepback, and higher wing load- 
ing have, in many oases, led to unsatisfactory dynamic lat- 
eral stability. Airplane designers are therefore finding it 
necessary to make such calculations in connection with the 
design and modification of airplanes. In many eases these 
calculations are difficult to perform for designers who have 
had no previous experience in theoretical stability work 
because most of the published theoretical analyses are not 
presented in a form that is especially suited to the compu- 
tation of dynamic stability. The estimation of the stability 
derivatives required in dynamic stability calculations has 
also been found to be difficult in many cases. Although 
theoretical and experimental data on these derivatives have 
appeared in numerous publications, no single publication has 
presented methods for estimating the derivatives for all 
types of airplanes. 

One approach to a presentation of methods of calculating 
stability and estimating stability derivatives in a form suit- 
able for use by designers was made by Zimmerman in 
reference 1. Although this report lias proved to be of valu- 
able assistance to designers in making dynamic stability 
calculations, recent trends in airplane design have caused its 

i Supersedes MAC A TN 2409, “Summary of 
and Marion 0. McKinney, 1951. 


usefulness to be seriously limited. For example, the equa- 
tions of reference 1 do not include the produet-of-inertia 
terms which have been shown by recent studies to be very 
important in some cases. (See references 2 and 3.) More- 
over, the calculation of the time histories of lateral motions, 
one type of calculation that has been the subject of increasing 
interest in the lust few years (references 4 to 7), is not covered 
in reference 1. The methods of estimating stability deriva- 
tives presented in reference 1 are also limited because they 
apply only to airplanes having unswept wings with an aspect 
ratio of 0 and operating at. speeds at which compressibility 
effects are negligible. The purpose of the present report is 
to extend the methods of reference 1 to include the methods 
of computation which are of current interest to designers 
and to include methods of estimating derivatives for con- 
figurations and flight conditions which are now being 
considered. 

This report summarizes and reduces to simple straight- 
forward steps methods for computing the time histories of 
lateral motions, the period and damping of these motions, 
and the lateral stability boundaries. Existing methods of 
estimating stability derivatives for a variety of airplane 
configurations are summarized and, in some cast's, simple 
new empirical formulas arc presented. Reference is also 
made to reports presenting experimental data that, should be 
useful in making estimates of these derivatives. 

SYMBOLS 

All forces and moments are referred to the stability 
system of axes which is defined in figure 1. The following 
definitions apply to the symbols except where they are 
otherwise defined: 

mass of airplane, slugs 
wing area, square feet 
wing mean chord, feet (i b/A ) 
wing span, feet 

span of that part of wing that has tip 
dihedral, feet 

tail length (distance from center of 
pressure of vertical tail to center of 
gravity, measured parallel to longi- 
tudinal stability axis; values of l 
must be calculated for each angle 
of attack), feet 

and for Estimating Lateral Stability Derivatives” by John J\ Campbell 
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1' i<*t ; KE 1. — Th« stability system of axes. Arrows intimate positive directions of moments, 
forces, ami angles. This system of axes is defined as an orthogonal system having the origin 
at the center of gravity and in which the Z-axis is in the plane of symmetry ami perpendicu- 
lar to the relative wind, the A'-axis is in the plane of symmetry and perpendicular to the 
Z-axis, and the V-uxis is perpendicular to the plane of symmetry. At a constant angle of 
attack, these axes are fixed in the airplane. 


A 

A 

X 


r 


IV 

i 

V 


kz 


0 


k x 


k z 


A'v, 


*-v„ 

b 

Kx » b 




kz 

b 



h 

w 


X 


d 


average fuselage height at wing root, 
feet 

average fuselage width at wing root, 
feet 

vertical distance of quarter chord of 
wing root chord from fuselage 
center line, positive downward, feet 
nondimensional time parameter based 
on span (1 7 tjh) 

longitudinal distance rearward from 
airplane (-enter of gravity to wing 
aerodynamic center, feet 
longitudinal distance from leading 
edge of vertical tail chord to hori- 
zontal tail aerodynamic center, feet 
(see fig. G) 

vertical distance from horizontal tail 
to base of vertical tail, feet (see 
fig. 6) 
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height of center of pressure of vertical 
tail above longitudinal stability axis; 
values of 2 must be calculated for 
each angle of attack, feet 
aspect ratio 

swccpback of wing quarter-chord line, 
degrees 

taper ratio (Tip chord/Root chord); 
also, differential operator in Laplace 
transform 

dihedral angle, degrees (see sketch of 
fig. 9) 

dihedral angle of wing tip, degrees 
time, seconds 
airspeed, feet per second 
radius of gyration about principal 
longitudinal axis of inertia, feet 
radius of gyration about principal 
normal axis of inertia, feet 
radius of gyration about .Y-axis, feet 

(V^ v ,r’ <,,,s ' * +i‘* 0 2 s ' n2 v) 

radius of gyration about Z-axis bud. 
(V^V* *? 'f £.v 0 2 sin 2 1 ?) 


product-of-inertia factor 
((kz*—k x *) sin rt cos v ) 


angle of attack of principal longitudi- 
nal axis of inertia, degrees (sec fig. 2) 
angle of climb, degrees (see fig. 2) 
angle of attack of longitudinal laxly 
axis, degrees (see fig. 2) 
angle between principal longitudinal 
axis of inertia and longitudinal body 
axis, degrees (see fig. 2) 
air density, slugs per cubic foot 
angle of bank, radians 
angle of yaw, radians 
angle of sideslip, radians 
rolling velocity, radians per second 
(d<t>/dt) 
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Figure 2.— System of u 

xos ;m<l angular relationship in flight. Arrows indifiatii positive 
direction of smiles, v =«■ t. 

r 

yawing velocity, radians per second 
((It! <lt) 

4> 0 

initial angle of hank, radians 

^ 0 

initial angle of yaw, radians 

fio 

initial angle of sideslip, radians 

(IHh 

nondimensional initial rolling velocity 

(// 0 /V/ (J ) 

(Ctft), 0 

nondimensional initial yawing velocity 

(</ \f/ 1 <1 <7 ) 

II 

Rout.lds discriminant or real part of 
complex root /? + /£ 

I 

imaginary part of (complex root R-\-Ii 

A,B,(,D,K 

coefficients of the (characteristic bi- 
quadratie e<| nation 

I\,I\ ...p. 

factors of the B y C, and I) coefficients 

^1 A2A3J V 

roots of characteristic biquadratic 
equation 

D 

differential operator ( djda ) 

p 

period of the lateral oscillation, seconds 

T m 

time to damp to one-half amplitude, 
seconds 

T 

time conversion factor (m/pSV) 

a 

nondimensional time factor (f/r) 

M 

relative density factor ( m/pSb ) 

L c 

impressed rolling moment, foot-pounds 

N c 

impressed yawing moment, foot- 
pounds 

Y c 

impressed lateral force, pounds 

c h 

impressed rolling-moment coefficient 

C„ c 

impressed yawing-moment coefficient 

Cr t 

impressed lateral-force coefficient 

Cl 

lift coefficient (Lift /qS) 

C D 

drag coefficient (DnigfqS) 

c, 

rolling-moment coefficient 
(Rolling moment /qSb) 

c n 

yawing-moment coefficient 
(Yawing moment /qSb) 

Cy 

late ral -f o rce coeffi c ien t 
(Lateral for ve/qS) 

g 

dynamic pressure, pounds per square 
foot Qp V s ) 


C,= 


dr. 
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Vr- 


<!rr 

4n 


l - ^ C ' c 
c 2K X 2 

^'n c 


Cr c 

Vc- 4> 



increment in C n produced by lift and 
induced-drag forces 


increment in C n produced by drag 
not associated with lift 

II 

horizontal tail 

a 0 

Subscripts: 

section lift curve slope 

vying 

wing 

fm 

fuselage 

tail 

used to designate vertical tail 


used to designate design under con- 
sideration 

data 

used to designate design for which 
force-test data are available 

exp 

experimental 

V -tail 

V-tail 

e 

effective 

II 

horizontal tail 


CALCULATION OF LATERAL STABILITY AND 
RESPONSE 

Various types of calculations may be performed to indicate 
in some way the stability of an airplane or the response to 
gust disturbances and control manipulations. The calcu- 
lations most commonly made are calculations of time his- 
tories of disturbed motions, period and damping of the free 
motions, and spiral and oscillatory stability boundaries (lines 
of neutral damping of the spiral mode and of the lateral 
oscillations). Step-by-step procedures for performing these 
types of calculations are explained in the text and derivations 
and additional pertinent material are presented in appendixes 
A to D. 

The period and damping calculations are the easiest of 
the three types to perform. For this reason, and because the 


dynamic lateral stability of airplanes is at present specified 
in the flying-qualities requirements in terms of the period 
and damping of the lateral oscillation, period and damping 
calculations arc probably the- most commonly performed. 

Recent dynamic stability work has indicated, however, 
that the period and damping characteristics of tin* free*, 
motions of an airplane a re not always a sufficient indication 
of whether the dynamic behavior of an airplane following 
various types ol disturbances will he considered satisfactory. 
For this reason the calculation of time histories of the 
motions ol airplanes is becoming more common despite 
the fact that these calculations are fairly laborious. The 
increasing use of automatic computing machines has also 
made the calculation of motions more popular. 

For many years, calculations of stability boundaries were 
the type oi calculation most commonly performed. In 
recent years, however, stability boundaries have not been 
considered to give an adequate indication of stability. 
Since boundaries are useful in some rases, however, (for 
example, for quirk approximation of the effects of changes 
in dihedral and tail area) the methods of calculating the 
spiral and oscillatory stability boundaries are described 
herein. Lines of constant period and damping of the lateral 
oscillation are related to stability boundaries (lines of 
neutral stability). In some cases these lines of constant 
period and damping may prove more useful than boundaries. 
Since no extensive use 4 has been made of lines of constant 
period and damping, however, tin* methods of calculating 
these lines (presented in references X and 9) are not given 
in the present report. 

The (‘(piations and methods of calculation presented in 
the present report deal specifically with the inherent motions 
of airplanes for the case of three degrees of freedom (roll, 
vavv, and sideslip) and linear stability derivatives. In order 
to perform similar calculations for cases involving additional 
degrees of freedom, nonlinear derivatives, or autopilots with 
lime lag, special equations are required. The methods and 
equations for treating these cases arc presented in references 
10 to IX. Additional degrees of freedom for the* ease of free 
controls are treated in references lfi to IX and for tin* case 
of fuel sloshing are treated in reference 10. The* use* of 
nonlinear derivatives in stability calculations is covered in 
reference 11. Methods of treating the effect of autopilots, 
including the effect, of time lag in the autopilot, art 1 presented 
in references 12 to 15 and 19. 

For some cases the effects of aerodynamic time lag art 1 
important. There are two different sources of such lag: 
(1) the time required for an aerodynamic impulse to travel 
from one component of the airplane 1 to another (for example, 
the time required for a change in sidewash at the wing to 
reach the tail a phenomenon commonly referred to as lag 
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of sidewasli) ; and (2) the time required ior the growth and 
decay of the aerodynamic loads on the airplane components. 

For both of these cases the time-lag effects usually become 
increasingly important as the period of the lateral oscillation 
decreases. The effects of the first type of time lag can be 
accounted for in some cases by modification of the stability 
derivatives. For example, the efleet of the lag of sidewasli 
on the derivative (\, T is discussed subsequently under the 
section on “Kstimation of Lateral Stability Derivatives. 

In many eases, however, both type's of time lag will require 
special stability equations. No general treatment- of these 
cases has been published but an indication of the method of 
treatment may be obtained from the treatments of autopilot 
lag in references 111 ami 15. 

CALCULATION OF PERIOD AND DAMPING 

As pointed out in references 1 and 2, the period and damp- 
ing of the various modes of the lateral motion may he calcu- 
lated from the roots of the characteristic equation 

A X 4 + ttX :i + C\ 2 -\ D X +K = 0 
by the equations 



and 

log, 2 0.603 

^ 1 / 2 = - ,, r « — H r 

where R represents a real root X or the real part of a complex 
root \ — II Ai Ii and / represents the imaginary part of a 
complex root. Negative values of 7 1/2 represent, the time 
required to double amplitude for unstable modes of the 
motion. 

The values of the coefficients A, H t D, and E may be 
obtained by the method given in steps 1, 2, and 3 of the sec- 
tion on “Calculation of Motions.” If the period and time 
to damp are to he calculated for a number of related cases, 
however, the values of the coefficients /l, R, (\ D, and K may 
be more conveniently calculated by a tabular procedure such 
as that shown as table I for making boundary calculations. 

Methods of determining the roots of the biquadratic 
characteristic equation are presented in appendix C. 

CALCULATION OF MOTIONS 

Calculation of the lateral motions of an airplane involves 
the integration of three simultaneous differential equations 
(see appendix A) to obtain a general solution in terms of the 
mass and aerodynamic parameters of the airplane. The 
general equations, once obtained, can then be used to obtain 
numerically the motions of any airplane in terms of the 
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variation with time of the angles of bank, yaw, and sideslip or 
some function of these angles such as rolling or yawing 
velocity. Various methods, such as those given in references 
20 to 22, are of course available for integrating the differential 
equations. Since the problems met in airplane dynamics are 
fairly complex, however, many of these methods are not 
suitable because of the difficulties of computation that arise. 
The method given in reference 4 (based on the Heaviside 
operational calculus) is satisfactory ior calculating tin* lorced 
motions following application of external forces or moments 
hut, without modification, this method cannot be used to 
calculate the motions resulting from initial displacements 
in bank, vaw, or sideslip or from initial values of rolling or 
yawing angular velocity. A solution based on the Laplace 
transformation is more satisfactory than that based on the 
Heaviside operational calculus because it permits direct 
calculation of the free motions following any initial condition, 
in addition to calculation of the forced motions following 
application of external forces and moments. The application 
of the Laplace transformation to the calculation of lateral 
motions is outlined in appendix B. The material presented 
in this appendix is similar to the work presented in references 
5 and 6 except that the mass and aerodynamic stability 
derivatives have been combined as shown in appendix A 
to reduce the number of arithmetical and algebraic, processes 
required in numerical solutions. 

The process of calculating the motions is presented as a 
scries of simple though lengthy arithmetical and algebiaic 
steps so that an understanding of the calculus involved in 
solving the differential equations is not required. 1 he 
method as shown is suitable for calculating the motions as 
variations of <f> 7 (3, />, and v with time for the case of the 

free motions following initial angular displacements (<t> o f 
\f/ 0y and j8 0 ) and angular velocities (D<t>)o and (7ty) 0 an *l f° r ^ 10 
case of the forced motions resulting from constant impressed 
forces and moments (L c , N e , and Y c ). These are the cases 
for which motions are usually calculated. It is also possible 
to calculate the motions resulting from impressed forces 
and moments which are arbitrary functions of time by the 
methods explained in references 0 and 7. 

MOTIONS RESULTING FROM INITIAL ANGULAR DISPLACEMENTS AND 
ANGULAR VELOCITIES AND FROM CONSTANT IMPRESSED 
FORCES AND MOMENTS 

The six steps involved in obtaining a specific solution for 
the lateral motions of an airplane are: 

Step 1: Determine values of the following parameters: 

(a) Mass characteristics 

m, A\v n , kz Q , 1 ?, and p 

(b) Geometric characteristics 


S and b 
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(c) Flight conditions 


r, C L} and 7 

(d) Aerodynamic stability derivatives 

C v (\- f , f\, C", and Cy t 

1 he methods of determining the values of the aerodynamic stability derivatives are given in subsequent sections of this report. 
In cases where impressed forces and moments are used as disturbances, determine the values of the factors 

<\> 

that are appropriate, to the particular problem. 

Step 2 : Prom the known {actors, evaluate the following parameters which are the stability derivatives in the form in 
which they are used in the calculation of motions: 


K - K ™ 

A, "Xv* 

K z ~ 

ft lr 
1 i 

il 

m 

~ p$b 

M 

<N 

II 

“"'2 K? C "* 

!E-\ <\. 


^=4AV^ 

, 1 r 

" 4AV "* 



l,= 4K x 2 Cl ' 

,>T ~4K~z* ( ”’ 

y,= 4p ( v ' 


when impressed forces and moments are 

used, (evaluate 



c 2Kx* ( l * 

v — r> 

c 2 K z 2 

?'•= ! r > . 



The values of Kx~, &z 2 y and K X z can be determined from the following expressions: 

K X *=K X * cos" rj-i K z * sin 2 rj 
K Z 2 —K Z{ ‘ 2 cos 2 ri +K X 2 sin 2 n 
&xz=(Kzq— K X( 2 ) sin t\ cos r) 


where 


k -*ia 

b 


Stop 3. Solve for the values of (lie appropriate ones of the following eocflicionls from equations ( 1 ) to ( 4 ): 
In all cases solve for the values of A, li, C, I), and E: 


where 


A=\-K l K, 

If— Ei — Ay s 

( — Ei Vs + E-, + 1’:, y P + P«y r — /*. 
E E s -tt+ E r , tan y + E- 


E— C.i -?r + E a — tan 7 


j 


Ei=-l v -nAKin v +K 2 l T 

E-i=l v n T —l T n v 


( 1 ) 
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— lull, — IrUj) 

l\= 4 „n ? —hn v 

Pi=Kin-a—h 

l\=K.J 9 -np 

I’ 7 = — I'-ilh + 1 \y p -t - 1 \}h— 1\ 


The quantities /', to P 7 are factors of the coefficients B, C, 1), and K which are combinations of terms that occur frequently 
in calculations of motions resulting from initial angular displacements and velocities and which are consequently grouped 
together for convenience. 

Calculate the values of a 0 , o u . . . a r> when solving for the angle of hank </> or the rolling velocity p: 

a x = 4 >J 1 + (I)<t>) 0 A 

a 2 — <f> , f T — ft o I ■ ) ( / )<t > )|»( — A Us - t - K J T n r) (I ) ^ 1 11 r I r) “h fc w c K\ 

«. 3 =0 n ( 1 \ -J- tan y-'rPi'j-foB:, - t y tan y- B&~ KJ ]- n^/^+iD^U- BiVr+ 

I\,-\ Ki V rlh — hih) — h(" r T >/j) + n r (K[ >/,i + /,) //<• Ci 

«4 = [<#> ll / > 4 -tAo/':<+(C ^)„/ > fi — \y tan y + lln $ —nfi!h + n r ij 3 )+n c Qtih—h—Byfi) — lhl > i 
ft s = (— Uitt+iiJfl) tan t 


Calculate the values of 6 0 , b u . . . b s when solving for angle of yaw i or the yawing velocity r: 
b „= tf'o.l 

b-,= t 0 O—pJ\—(I)<t>UKJ,,—n p )+{ItyU-Ay»+Ki n i’--l,,)—lcK°.+nc 
bz=-<t> n I\ £+* (l\ y+ /'*- ) - f KJ p y 9 —n p y „) -i - 

Kplpl/jj + lp 1 /$) + /c(^2'/d+ p) — + ?//») — ?/c Cl 

b i ^[-<j> n I\+^ a P 3 —(D<l>)J\+(Di)J\] ( -Y+lc('n^yp—n,,y$)+n c (l p y 9 —hy p )—!/cP 4 

b;—(J c i‘n—nJn) ( ~y 


Calculate the values of c„, c t , . . . c 4 when solving for the angle of sideslip fi: 
c o” &o A 

Ci = <j>QA ~~ Lj r - 4 ^ tan 7 + / 3 oPi+(/-^ 0 )o^ 4 .?/p (^ < / / )o^.( 2 /r 

C.., = 4> 0 Pi % + M —■ tan 7 + /J 0 C! + (/}<#>)() ^ — i^jd/r + ^^ + TVl/r - + — »r)jfpJ + 

(C^)o [yl^tan 7+Ada p y r -Adn,, -(A' 1 rG-fd!/,,]+«-^.Vr+AC+;/,>)+^(// 4 -l-/A ? / ; ,)+;/c/ J i 

c,= ^oC 2 Y + 'I'J’-I ~.y tan 7 + (D<j>)t, ( ~K 2 l p ~ tan 7 + % ^ tan y+K 2 l T y-n r ^r)+ 

(Dtx(K,n v ^y tan 7 -/ p ^ tan y-Kpn, ^+lr^£}+h(n v y r -n p -n T y p + ^- K ^ tanT ) + 

« c ^ — l P yrPlpplr y p ^1 tan 'Y^'h'f/eCo 

c i == l p (n p ^ tan 7 -»r ^r)+». (jr tan 7) 

234238— 53 2 
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Step 4: Solve for the roots X 1? X 2? X 3 , and X 4 of the biquad- 
ratic equation 

A\*+B\*+C\*+n\ + E =0 (5) 

where the values of the coefficients A , B, , . . etc. were 
gi vcn by the solution of equations (I). Methods of deter- 
mining t lie roots of the biquadratic equation are given in 
appendix C. 

Step 5: Use the coefficients obtained from equations (1) 
to (4) and the roots of equation (5) to solve for the following 
coefficients: 

Calculate the values of the factors A h A*. . . . A 6 when 
solving for the angle of bank <£ or the rolling velocity p: 


1 ”6ylX 1 5 +5/?X 1 4 + 4CX I 3 + 3/>Xrd-2/tX I 

_ ^0^2° 4“ d\ X2 4 ”b ®2X2'M“ 3 X 2^ “ U (l 4 \ 2 -) (l 5 
2 _ 6 3 X 2 5 + 5 # X 4 + 4 CX 2 3 4* 3 /? X 4 + 2 E X 2 


&oX3 5 ~b &i X 3 4 4~ fl 2 Xa 3 4' O3X3 2 4~ a 4 X:j -f- 

3 ~ 6 ,4 X^-jFsif X 3 4 + 4 (?X 3 3 + 3 X 3 “ 4“ 2 E X 3 

^0X4° + a 1 X 4 4 4~ g 2 X 4 3 4~ g 3 X 4 2 + a 4 X 4 H- a 5 

4 6^4 X 4 5 + 5 Z?X 4 4 + 4 CX 4 3 4~ 3 D X 4 2 4- 2 E X 4 


A h = 


a- 

E 


A 


E 


(«4-a 5 g) 


( 6 ) 


Calculate the values of the factors B h B 2 , . . . B & when 
solving for the angle of yaw \p or the yawing velocity r: 


» = 6 pXi*+ 6 iX 1 4 + 6 J XiH 6 aXi*+ 6 4 Xi + 65 "1 
6 M x4 -f- 5/^X4 4“ 4CXi 3 + 3/?X 1 2 -|- 2 A’X 1 

^oX 2 5 +6 i X 2 4 + 6 2 X4 + W+6 4 X*- 4 6, 

2 (j/1 X a 5 + 5 »X/+4<7X? 3 D\/ + 2E\> 

U __ 6 oX 3 5 + ^iX 3 4 + & 2 V 4 - 6 3 X 3 2 4 ^ 4 X 3 + 6 5 
6 >4 X 3 “' 4~ 5B \ 3 4 ~b 4 CX 3 3 4 3/)X 3 j 4“2 E X 3 

jj ^oX 4 5 + &iX 4 4 + 6 2 X 4 3 -b ^ 3 X 4 2 4 ~ 64X4 ! ^5 

4 ~ 6,4 X 4 5 + f> «X?'+4(7X 4 3 4-30 X 4 2 + 2 ^X 4 


*(*• ‘") 




(7) 


Calculate the values of the factors (\, C 2 , . . . C 5 when 


solving for the angle of sid(‘sli|) fi : 


CoX,” b c ! X ] 4 - 1 r>Xi 3 -l c 3 X] 2 4 c 4 X| 1 

1 6,1X4 ^ 5/yx, 4 4 ■4f ri X 1 3 + 30X 1 2 ' 4-2A T X, 

(i c iy \ 2 * 4- c 1 X4 4- c-X./ + c 3 X 2 2 4 c 4 X 2 

2 6 A X/' ! 5 /fX , 4 + 4 r 'X 2 :i -4 3 1) X 2 2 + 2 E X 2 

CoXj^ + CjXa 1 -!- CaXj/q ^3X3"-! c 4 X 3 
6 ,4 X 3 ' 1 i5 /iX 3 4 4 4 1 ’X3 3 4 ■ 3 1) X 3 " + 2 E\ 3 

^oX 4 5 4 _ CiX 4 4 d- c 2 X 4 * 1 -b<‘:iX 4 J 4- c 4 X 4 

4 «-4X 4 s - 5tfX 4 4 + 4C\ 4 M 30X/ + 2/i’X 4 


<V 


E 


If equation (5) has conjugate complex roots, the values of 
the coefficients (equations (6) to (8)) corresponding to these 
roots will be conjugate complex. In order to facilitate 
treatment of this ease it is convenient to establish some 
special notation. This special notation is explained in 
appendix 1). 

Step 6: The equations of motion are written in different 
form depending upon the roots of equation (5). If the 
characteristic equation has four real roots X 1? X 2 , X a , and X 4 , 
the general form of the equations of motion is used as follows: 


^=-4^1 ■b-4 a ***4--4,«'M A 4 e 9X *+A B tr + A 9 
1 p=B x v" K i ; B 2 e^-X-B^ j-B,e a ^ + B,<r + B 0 
0 - - ( V* K * + (\e' k * + (\e 9k * + (\ 

P — ( * 1 1 X 1 e ff x 1 , 4 2 X 2 (- 4~ A 3 X 3C ffX; i - f , 1 4 X ^ K * j .1 5) 

T ' I 


0>) 


r - 1 ( B x \^ 1 + B 2 \ 2 e OK * + B,\,e^ + B A \ 4 e' K * + B h ) 


If ? as is gtuuually the case, equation (5) has two complex 
roots and two real roots ft— II, X 3 , and X 4 ), the equa- 

tions of motion may be expressed as 


<P = K A <° {i cos (<r 1 4 - w 4 ) 4* + A \e° x * + A h <r 4- ,4 6 > 

\p — K h E r (*os (<t 1 4" B^e aK:i 4" B A e vx ^ I B h <j -f- B§ 

P cos Ul + uJ+f 4- (\e^ .\-C 5 


P—* T \ ft 2 4 Ee aR cos(a/ + w 4 4-tan~ l y^ + 

* 1 3X3^^^ 4~ 44X4^4- /ls^ 

T ~ t ^ ^ ( ° K <>os ^°’^"I" C0 *4-tan“ 1 ^>^) + 

/f 4 X 4 U x 4 + /^"| 


f 00) 


j 
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f'ALCt’LATIMi LATEHAL STABILITY A XI) ESTIMATIN' (} STABILITY DE11I VATIVES 


wllCTP 


2 ;/,y- 

t he 

W,t“ 

= tan 

2 \ /in" 

I/,, 2 

OJ/i 

-tan 

2\7i’ r -- 

\ hr 

0i c 

tan - 


, rr 

Ri 

, Ih 
R„ 

h 
lie. 


(10m) 


ami and /, arc defined in appendix 1). 

If there are four complex roots (I! ( li, R — h, R' T'i, and 
R ' — /'/), the equations are 

(f} — K A e' ,H cos (<7 1 e" 1 *' cos (al' ' 

A,- t a +a1 6 

}p=K„e ,Tl1 cos (a / d- co/t) - K n’ < aR ' cos (cr /'+««')- f- 
Ito -f- it - 

P=K C t‘' tt cos (cl + Uc )+K c ’f?* cos (cT +«</) + <?* 


[/v , \R i \ I 2 e° h cos 

[<rl ~ o>a -f tan 1 -(- A 5 + 

K A \R ri 1 /'V"' cos ( 

<r/'4 CO,/ ft an -1 J 

,■ ' | /V"cos 

(^(rl 03 h~\~ tan 1 | 

l n r H ‘ cos ( 

* /' + «,/ + tail' 1 

when' 


AY 2 \ R,i 2 r /./ 2 

o}a — tan 1 

AY=2 \ AY 2 - Ib 2 

j f 

t a _i j n 
avZ—tan 1 , 

J c f 

AY 2^R c n \ I7 2 

u>c — tan 1 ft f 


( 11 ) 


Hi la) 


The coefficients K A , K„, K c , o>a, w,„ and co c are defined in equa- 
tions (Ida) and R A , I A , R/, and /./ are defined in appendix D. 

Solve the appropriate ones of these equations of motion 
(equations (!)), (10), or (11)) by substituting values of tins 
nondimensional time factor tr in the equations and solving 
for <t>, f , P, p, or r. 

MOTIONS RESULTING FROM ARBITRARY DISTURBANCES 

The motions resulting from arbitrary forcing functions 
can bo obtained from the motions resulting from constant 
impressed forces and moments by the methods explained in 
references 6 and 7. 

A very useful method of obtaining the motion resulting 


from various abrupt gust and control disturbance's is given 
by Jones in reference 7. In this report it is pointed out that, 
although the component motions of an airplane must be 
calculated simultaneously (that is, by simultaneous difler- 
ential equations), the effects of component disturbances may 
by the principle of superposition be calculated separately 
and later added in any desired proportion. Thus, if a given 
rolling moment causes a 20° bunk in 1 second and i! a given 
yawing moment causes a 5° bank in 1 second, the combined 
effect of both acting simultaneously will be a 2f)° bank in 1 
second. Jones also points out a somewhat similar fact with 
regard to the effects of disturbances that are not applied 
simultaneously. This fact is that, il a given disturbance 
which arises at the time t 0 is later augmented, the effect 
of tin* increment of disturbance will run its course inde- 
pendently of the effect of the original disturbance, hor 
example, in a problem involving the correction for a gust 
disturbance by a manipulation of the control, the motion 
produced by the gust disturbance can be calculated inde- 
pendently and the motion caused by the assumed corrective 
control manipulation can he added to it at any desired point. 
This example is illustrated graphically in figure 3. 

The principle of superposition may be applied analytically 
as well as graphically. The analytical application which 
makes use of Carson s integral or Duhamol’a integral is 
described in references 7 and 23. This method is useful for 
calculating the motions resulting from impressed forces and 
moments which are arbitrary functions of time. By applica- 
tion of those methods, the solutions for constant impressed 
forces and moments can be used to obtain new solutions for 
any arbitrary variation of impressed forces and moments 
with time which can be expressed by a mathematical formula. 
Some simple variations of impressed forces and moments 
with time and their Laplace transforms are given in reference 
B. The transforms for any other function for which trans- 
forms have been worked out may be found in tables of 
I ja p 1 ace tr ai l s form s . 

CALCULATION OF STABILITY BOUNDARIES 

OSCILLATORY STABILITY BOUNDARIES 

As pointed out in the preceding section of this report, the 
degree of stability of the uncontrolled mot ions of an airplane 
is indicated by roots of the characteristic equation 

A X 1 I- bV + + D\ + E= 0 

For stability the real roots or the real part of the complex 
roots of the characteristic equation must be negative. A 
useful discriminant for determining some of the character- 
istics of the roots in stability work is Kouth’s discriminant 
ft (ft=BCD—AJ) 2 —B 2 E). The use of this discriminant in 
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dynamic stability analyses has been point ad out in many 
reports, for example, references 1, 2, ">, 21. and 24. Kouth 

has shown (reference 20) that, if It and (lie coefficient E are 
finite, the necessary and sufficient conditions that the real 
roots and the real parts of the complex roots should be 
negative arc that every coefficient of the biquadratic and 
also It should have the same sign. Kouth also showed that 
when //=() and /> and 1) have the same sign there are a pair 
ot complex roots with tin* real parts zero. Since the value of 
the real part of a complex root indicates the stability of an 
oscillatory mode of the motion of an airplane, the lateral 
oscillation is neutrally stable when It ^ 0 and the coefficients 
7> and I) have the same sign. Oscillatory stability boundaries 
can be determined, therefore, by solving the equation It— 0 
and checking to determine whether the signs of /> and 1) are 
the same. 






(a) Gust disturbance. 

(b) Control manipulation. 

(c) Gust disturbance and corrective control manipulation. 

Fzucrk 3.— Illustration of superposition of motions t<> determine effect of arbitrary 
disturbances. 


Since two ot the most important stability derivatives af- 
fecting lateral stability are the directional stability derivat ive 
* UU( I I he (41cct ive dihedral derivative (\ y houndaries for 
muitral oscilhilory stability are usually calculated as a func- 
tion of these two derivatives as illustrated in figure 4. These 
calculations arc generally curried out by tin* method shown 
111 table 1. This table contains a numerical example* and 
step-by-sl ep instructions for using the table. Tin* results of 
this numerical example an* plotted in figure 4. The pro- 
cedure illustrated in table 1 is first to assume values of tin* 
independent variable ( \ f(j to cover the range* for which the 
boundary is required. I'lie values of all the other mass and 
aerodynamic stability derivatives except (\. an* llien esti- 
mated. The value of ( is generally assumed to have* been 
varied by varying the size of the vertical tail and, con- 
sequently, tin* tail contribution to each of the other stability 
derivatives varies as ( is varied. Tin* values of the coef- 
ficients A, />, (\ l) r and E and then It are calculated as 
fund ions of / rf : 


The value's ot f ;i corresponding to the assumed values of (\ 
for the condition of neutral oscillatory stability arc next 
obtained by solving the expression II 0 which is a quadratic, 
in I ;i t bat is of the* form 

11 \! / + Vi I d v'] ■ 0 

Mindly, tin* values ot ( corresponding to tin* assumed 
values of ( \ ffj an* obtained from tin* values of /#. 

1 he value's of f ;i w hich satisfy t he expression It 0 must In* 
ch(*ck(*d to determine whether they satisfy tin* other con- 
dition for neutral oscillatory stability that the* sign of the 
coefficients /> and I) must hi* the* same*. This check can be* 
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pi'ii'oriiii'd readily by substituting the values of which 
satisfy 7/=0 into the expression for I) which is a linear 
cquat ion of the form 

I )-- ?/o 

Thus, the sign of I> is determined. The sign of /> is a con- 
stant for any given value of (\ iQ and is almost invariably 
positive sinct* tin* three predominant terms of /> contain the 
derivatives (\ , (\, , and (\- which in all practical cases 

p r d 

contribute a positive increment to the value of />. 

Since two values of (\^ satisfy the condition II— 0 for each 
valu(‘ of (' n , the // 0 curve lias two branches. As pointed 

out in reference 24, one of the branches ot the /T O curve 
generally represents an oscillatory stability boundary and 
the other branch represents a line of numerically equal real 
roots with opposite signs. (See fig. 4.) if neither of the 
values of (\ which satisfy the expression II — 0 for a par- 
ticular value of (J n is found to represent- a point of neutral 
oscillatory stability, tin* lateral motion has no oscillatory 
mode for that value of . If both of the values of Ci fi 
which satisfy the expression. /T O an* found to represent, 
points of neutral oscillatory stability, the lateral motion 
lias two oscillatory moth's. In this case, since tin* boundary 
D—i) represents the line of infinite period, the branch of the 
}>~ 0 boundary which lies close to tlu* //--() boundary is 
usually tilt' boundary for neutral stability of the longer period 
of the two oscillatory inodes. A detailed discussion of the 
significance of the stability boundaries and the regions formed 
by these boundaries is given in reference 24. 

In calculating stability boundaries for a specific airplane 
a complete solution such as that explained in the preceding 
paragraphs should be made. For general studies of stability, 
however, approximate oscillatory stability boundaries may 
lx* calculated much more simply by the methods shown in 

reference 24. 

As pointed out previously, methods of calculating lines of 
constant period and damping of the lateral oscillation are 
presented in references S and 9. 

SPIRAL STABILITY BOUNDARIES 

Spiral stability boundaries, like oscillatory stability bound- 
aries, are usually determined as a function of the directional 
stability derivative' fk and the effective dihedral derivative 

•• p 

(\ as illustrated in figure' 4. As pointed out in reference 1 ; 
neutral spiral stability occurs when the K coefficient of the 
characteristic equation is zero (/'/==()). A spiral stability 
boundary can be easily obtained from this relation. If 
expressions for K (in terms of In) corresponding to several 
values of C n have already been obtained in the process of 
calculating an oscillatory stability boundary, the equations 
formed by setting these expressions for K equal to zero can 
b(' solved for the values of In (and hence (\ p ) corresponding 
to the assumed values of . If the values of K have not 
already been obtained in the process of calculating an oscil- 
latory stability boundary, a spiral stability boundary for the 
level-flight, condition ( 7 ^ 0 ) can he calculated simply from 


the ('(illation 


Values of (\ l& are assumed within the range for which the 
boundary is required. The values of (\ r and (\ r corre- 
sponding to each value of (\ % an* then determined. The 
tail contributions to these derivative's generally vary with 
(\ since (\ e is usually a-ssumed to la* yaried by changing 
the size of the vertical tail. 

ESTIMATION OF LATERAL STABILITY DERIVATIVES 
GENERAL REMARKS 

Methods of estimat mg tlx* lateral stability derivat ivos have 
been presented in numerous publications but no single report 
has contained information for estimating the contribution ot 
all principal airplane components to all the derivatives for 
airplanes having any sweep angle or aspect ratio. In the 
present report, an approach to such a presentation is made 
by the coordination of and reference to existing estimation 
methods, hv reference to publications containing data which 
should be useful in making estimates, and by the suggestion 
in some cases of simple' new empirical formulas. Detailed 
estimation methods arc presented for low-subsonic-speed 
conditions but only a brief discussion and a list of references 
are given for transonic- and supersonic-speed conditions. 
In general, the estimation methods presented should be ex- 
pected to yield only fairly accurate values suitable for making 
first approximations of dynamic stability. This limitation 
applies especially to the cases in which tlx* derivatives are 
based completely on theory. 

For convenience, the references that should he useful in 
estimating the stability derivatives are presented in table II. 
The references are grouped according to the speed range 
covered (subsonic or supersonic) and according to the deriva- 
tives presented in each report. The references for the sub- 
sonic case (references 1 and 25 to 97) are further divided into 
two groups one including reports which contain estimation 
methods and the other including reports which contain 
experimental data that should bo useful in making estimates 
of derivatives. The references for the supersonic case (refer- 
ences 98 to 1 18) are subdivided according to wing plan form. 

The following sections covering the estimation of the nine 
stability derivatives are divided into three groups according 
to the type of derivative -sideslip derivative's (/ * t 8 ) } 
rolling derivatives ( r p h «d <1 yawing derivatives 

(( r n r} (\ r r y> ). The derivatives C Yp and (\ T have usually 
been neglected in making dynamic lateral stability calcula- 
tions because theory indicated that for unswopt wings C Yp 
and (\ r were zero. Recent experimental data, however, 
have indicated that both swept and unswept wings produce 
measurable values of these derivatives (references 25, GO, 
and 89). Since the vertical tail contributes to ( y y p and (\ Ti 
it appears desirable to estimate these derivative's and to use 
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them in tin* calculations of stability unless it is ost ablishctl 
that for the case in question the effects of (' y and (\ r on 
stability are negligible. For these two derivatives, only the 
effect of the wing and vertical tail need to be considered. 

The methods of estimating the rolling and yawing deriva- 
tives presented herein were obtained from theoretical treat- 
ments based on the assumption of steady rolling and yawing 
and from experimental data obtained principally from tests 
made under conditions of steady rolling and yawing. The 
only information that applies directly to the oscillatory 
case is a limited amount of data, on (\ f obtained by oscillation 
techniques. When calculations are made in which the 
oscillatory mode is the subject of interest, some consideration 
should he given to correcting the derivatives based on steady 
rolling or yawing to account for differences in the derivatives 
that are likely to exist as a result of differences between 
the oscillatory motion and the steady rolling and yawing 
motion. For example, tin* data of reference 85 have indi- 
cated that, for flap-extended or power-on conditions, fairly 
large differences might exist bet vveen the values of the tail 
contribution to (\ r for the steady yawing and yawing oscil- 
lation (‘uses. At present, little information is available for 
correcting the values of (\ r for the steady yawing ease to 
apply to the 1 oscillatory ease and, unfortunately, little or 
no information is available for correcting the other stability 
derivatives. 

Since most wind-tunnel force-test data that are likely to 
be used in making estimates of the stability derivatives 
are probably for much lower Reynolds numbers than those 
for the full-scale airplane, some adjustments to the data 
are usually required to account for the differences in Reynolds 
number. The effects of Reynolds number should be con- 
sidered in the cases of all the derivative's, especially those 
which arc 4 estimated by methods that involve the use of 
force-test data. Methods of correcting for Reynolds num- 
ber effects for some of the derivative's an 4 discussed in the 
following sections which cover the estimation procedures. 
In the 4 cases where* the Reynolds number effects are not 
discussed, it ran be 4 assumed that any abrupt variation in 
the derivatives near the stall for low-scale data will also 
be present for the full-scale airplane but will probably occur 
at a higher lift coefficient because 4 of the higher maximum 
lift coefficient, of the airplane. An indication of the lift- 
coefficient range over which the theory may not he expected 
to give reliable values of stability derivative's for the* full- 
scab 4 airplane can be obtained from large-scale drag data. 
The analysis of reference 4 89 indicates that the variation 
of the derivative's with lift eocfficicnt is different from the 
theoretical variation at lift coefficients above that at which 
the drag elue to lift increase's abruptly from the ideal value 

<VhA. 

The effects of Mach number and power are not treated 
in the sections on the individual derivative's hut are dis- 
cussed briefly in separate sections. A detailed treatment 
of these effects, including design formulas and charts, was 
considered bevonel the scope of this report . 


THE SIDESLIP DERIVATIV ES C v , C , C, 

1 a n a l a 

No satisfactory purely theoretical methods have 4 yet been 
developed for obtaining accurate estimates of tin 4 sideslip 
derivatives and for a complete airplane, pri- 

marily because of large interference effects between the 
various airplane components and because of large, and often 
unpredictable, variations of 1 he derivatives with angle of 
attack, t ort unat oly, these derivatives can be obtained 
from conventional wind-tunnel force-test data. Such experi- 
mental data are essential to the accurate determination of 
sideslip derivatives. It is, of course, highly desirable to 
have force-test data for the exact- airplane design under 
consideration, but reasonably accurate estimates can usually 
be made by correcting the force-test data for a generally 
similar design. The methods of correcting the force-test 
data on a similar design for use in the (‘use under consideration 
are covered in tin 4 following sections. In the formulas 
presented, the subscript word “design” is used to designate 
the design under consideration and the subscript word 
“data” is used to designate the similar design for which 
force.- test data an* available. 

Force-test data should be used to determine the effect on 
tin 4 sideslip derivative's of such airplane components as 
leading-edge high-lift devices, stall-control devices, trading- 
edge flaps, nacelles, external stores, canopies, and dorsal and 
ventral fins. The effect, of leading-edgi' high-lift devices is 
usually merely to extend to a higher lift coefficient the same 
variation of the derivative with lift coefficient as for the plain 
wing. Trading-edge (laps often have large effects on the, 
contributions of both the wing and the vertical tail to the 
sideslip derivatives (references 40 and 71); and since these 
effects an 4 not easily estimated, it appears that in these eases 
use of force-test data is essential. The addition of nacelles 
and external stores generally lias been found to decrease the 
directional stability factor (\ ti slightly. The results of a 
limited amount of research to determine the effect, on the 
sideslip derivatives of the size and shape of canopies has 
been reported in references 49 and 75 hut these results are 
inadequate for making accurate predictions of the effects 
of canopies. The effects on the sideslip derivatives of dorsal 
and ventral fins are usually small at the small and moderate 
angles of yaw that are generally considered in stability 
calculations. (See references 48 and 73.) 

r> . 

In estimates of the lateral force due to sideslip derivative 
force-test data for the design under consideration 
should be used whenever possible. If such data arc not 
available, data for a similar design can be used and corrected 
as follows: 

Wing-fuselage. Since the wing-fuselage contribution to 
C Yfi is usually relatively small compared with that of the 
vertical tail, great accuracy is not required in estimating this 
factor. This contribution may be estimated as follows: 

(1) Wing: If tlie wings of the two designs are generally 
similar, the difference in (\^ can he considered negligible 



METHODS 


FOR CALCULATING LATERAL STABILITY AND ESTIMATING STABILITY DERIVATIVES 


13 


and no correction is necessary. I lie theory of reference 2o 
docs not appear to he suitable for use in estimating * 

(2) Fuselage: If the two fuselages are similar in shape, 
the difference in C Y(j can probably be estimated satisfac- 
torily by correcting for the difference in the relative size of 
the fuselage and wing for tin 1 two airplanes. It appears, 
however, from table X of reference 71 unlikely that a reliable 
prediction of C Y can be made directly from the geometry 

1 ... s-, 

of the fuselage. Some additional data on ( Y are pre- 
sent(‘d in reference 79. Experimental data from other 
investigations have shown that differences in fuselage 11 cross 
section can cause* very large* differences in the variation of 
C Y with angle of attack. For example, in the ease of a 
flat fuselage with the major cross-sectional axis horizontal, 
the sign of (* y has been found to reverse at moderate and 
high angles of attack. lorce-test data are essential for 
making estimates in such cases. 

(:V) Wing-fuselage* interference*: For low-wing or high-wing 
configurations, wing-fuselage interference causes the value 
of C Y to be greater than that obtained by adding the con- 
tributions of tin* wing and fuselage. (See references 37 and 
40.) If the vertical location of the wing on the fuselage is 
generally similar for the two designs, however, any correction 
for a difference in this interference factor can he neglected. 

Vertical tail. Accurate estimates of C Y are necessary 
because this factor is used to estimate the tail contribution 
to several other derivatives. This factor is especially im- 
portant at low angles of attack because in this case the tail 
contribution is often much greater than the wing-fuselage 
contribution to all derivatives except (\ p . For this reason 
it is highly desirable to have tail-off and tail-on force-test 
data for the design under consideration or for a very similar 
design. Corrections to the data for a similar design can be 
made as follows: 

(1) Correction for differences in wing area, tail area, and 
tail lift-curve slope can be made by the following formula: 


( ( L & tali') . V 

V a tai l /design ^ O < 


* v- , \ / / ■> \ \ *tn>l /design 

\ ^ tai l / design \ & tail /data { ( £ S 

V a lail / data 


1 des ign 


( 13 ) 


The value of < can be obtained from figures 5 and 6 

a lail 

which are based on the theory of reference 34 and on the 
theory and data of references 28 and 35, respectively. The 
chart of figure 6 can he used to estimate the change in the 
effective aspect ratio of the vertical tail caused by the end- 
plate effect of the horizontal tail. It should he emphasized 
that for the best accuracy the charts in figures 5 and 6 should 
he used in conjunction with formula (13) for correcting 
existing force-test data and not for making a direct estimate 
of V Y . 

Ptail 

(2) In the ease of V-tails, the correction for C Yb , can be 

^ 1 M tall 

made as follows: 


( 1 f 




( * ) I ( 


(KC, S v .i„u sin 2 F) 


K ^d/ita 
k 'design 


where the terms (\ , T, and K are the same as given in 

reference 30 and arc defined as follows: 

(\ slope of the tail lift curve in pitch measured in the 
plane normal to the chord plane of each tail panel 
T dihedral angle of tail surface measured from A ) -plane 
of the tail to each tail panel, degrees 
K ratio of sum of lifts obtained by equal and opposite 
changes in angle of attack of two semispans of tail 
to lifts obtained by an equal change in angle of 
attach for the complete tail 

Values of the term K , which are usually about 0.7, can be 
obtained from reference 30. 

(d) Since large* differences in sidewash and dynamic pres- 
sure at the* tail can be caused by differences in wing plan 
form and wing location, use of experimental data for the 
specific design or at least for a design which has a closel\ 
similar wing-fuselage* combination and vertical tail location 
is extremelv desirable. No methods are* available* which 
permit accurate predictions e>f sidewash at the* tail, but the 
experimental elata of references 40, 50, and 71 can he* used to 
obtain some indication of the* variation in sidewash with 
vertical location e>f an unswept wing on a fuselage and the 
experimental elata of references 36 and 79 provide additional 
information on sidewash at the tail. Other experimental 
data indicate that the sidewash fields produced by highly 
swept, low-aspect-ratio wings or by fuselages of flat cross 
section can sometimes be strong enough at. high angle's of 
attack to reverse the effectiveness of a conventionally located 
vertical tail surface. Until a reliable method is developed 
for pre>dict.ing these large siele*wash effects, force-lost data 
appear to be the only means by which satisfactory estimates 
of C Ya can be obtained. 

Ptall 

C \ 

Although attempts have been made to develop methods 
for estimating the yawing moment due* to sideslip (static 
directional stability) derivative C n& (for example, references 
70 and 71), no reliable* method has yet been obtained. I he 
use of force-test data therefore* seems imperative. 

Force-test data for the design under consideration should 
be useel if available. If such data are not available, use* data 
for a similar design and correct as explainer! in the sections to 
follow. 

Wing-fuselage —The corrections for the wing-fuselage 
contributions are: 

(1) Correction for wing: From figure 7 (taken from refer- 
ence 25) the values of (C nfi /C L 2 ) wing for the design under 
consideration and for the design for which test data are 
available can be determined. The effect of differences in 
taper ratio can be neglected. (See references 61 and 68.) 
The difference between these values of C^/C^ 2 should then 
be added (with proper regard for sign) to the experimental 
data for the complete model. 

(2) Correction for fuselage: The formula 

/Fuselage volumeX / h \ 

V Sb A w) 


( 14 ) 


(15) 



* e tail (Hon) 
e taH ( h off) 
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for Hit* cast* of subsonic incompressible How. «=() '' Taken from reference 35. 


but should only be used as indicated to determine a correc- 
tion for force-test data. This correction method should 
not be used in the cases of high angles of attack when there 
are large differences in fuselage configuration. Force-test 
data are essential in sueli eases. 

(3) Correction for vertical location of the wing: If the 
designs arc generally similar, the correction for the vertical 
location of the wing on the fuselage can be neglected. (See 
references 37 and 40.) 

(4) Correction for eentor-of-gravity position: If the centcr- 
of-gravity position for the design under consideration is 
appreciably different from that for the design for which 
force-test data are available, the value of C n(t for the wing- 
fuselage combination can be corrected by multiplying the 
value of Cr 3 for the wing-fuselage combination by the 
distance between center-of-gravity positions (expressed in 


wingspans). 

Vertical tail. — Corrections to C n » for differences in 

PfOll 

Cy and tail length l/b can be made by the following 

Ptail 

formula: 


(<?«, 


) desi an 



( 10 ) 


The contribution of wing-tip fins to C nfi is treated in 
references 07, 72, and 87. 

In estimates of the rolling moment due to sideslip (effective 
dihedral) f force-test data for the design under considera- 
tion should be used. If such data are not available, data 
for a similar design can be used and corrected by the methods 
that follow. 



nn 
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Wing-fuselage. — The corrections for wing-fuselage con- 
tributions are: 

(1) Correction for wing: From figure 8 (based on reference 
25) the theoretical values of (\J(\ for the design under con- 
sideration and for the design for which data are available can 
be determined. The difference between these two theoretical 
values can then he added (with proper regard for sign) to the 
experimental data. Consideration should be given to scale 
effect, airfoil section, and surface roughness on the value of 

for highly swept wings. The lift coefficient- at which the 
experimental variation of C lff with lift coefficient departs 
from theory is greatest at high Reynolds numbers and for 
smooth wings with round leading edges. For wings with 
rough surfaces or sharp leading edges the effects of Reynolds 
number on art* usually small and low-scale wind-tunnel 
data can be used. For airplane's having very smooth swept- 
back wings with rounded leading edges, however, some cor- 
rection should he made for scale effect, when estimations are 
made from low-seal <5 wind-tunnel data. Since no rational 
method lias been developed for making such corrections, it is 
suggested that, for lift coefficients higher than that at which 
the experimental data departs from the theory, an average 
of the theoretical and low-scale experimental values he used. 
Conservative dynamic stability results will usually be ob- 
tained if the uncorrected theoretical values of (\ are used 
because these values are ordinarily greater (more negative) 
than measured values and because tin* larger negative values 
of <7*0 usually tend to decrease the* dynamic lateral stability. 

(2) Correction for wing dihedral: The effect of dihedral 
on is treated in references 29, 40, 52, 59, 08, and 81. 
Correction for the difference in dihedral between the two 
designs can be made by multiplying the incremental geo- 
metric dihedral angle (in degrees) by the factor f V obtained 

from figure 9. A plot of (', against aspect ratio for taper 

ratios of 1.0, 0.5, and 0.25 (obtained from references 59 and 
08) and a formula from reference 59 for correcting for sweep 
are presented in the upper portion of figure 9. The lower 
chart and formula in figure 9 (developed from reference 08) 
should he used in addition to the upper chart and formula 
of figure 9 to estimate the values of for the ease of a w ing 

with partial-span dihedral. Although this chart and formula 
apply directly only to wings with one dihedral break, they 
can be used to estimate the (\ for wings with two or more 

dihedral breaks by the method described in reference 08. 
Iho effect, of drooped wing tips and of wing-tip end-plates 
011 ^ {{S win a s ^ ou ^ determined by experimental data since 

no reliable estimation procedure for these effects is available. 
(See reference 07.) 

(d) Correction for wing-fuselage interference: Although tilt* 
contribution of the fuselage alone to (' t is usually negligible, 
the interference between the whig and fuselage can greatly 
alter the value of of the wing. This interference is such 
that a high location of the wing on the fuselage gives more 
positive effective dihedral (higher — C, ) and a hnv wing 
location gives less positive dihedral than a midwing position. 


This effect is treated theoretically in reference 69 and lias 
been studied experimentally in references 37 and 39 to 43. 
The following simplified expression for estimating the incre- 
ment in (\^ caused by wing-fuselage interference has been 

developed from the relationships presented in reference 69 
and in other sources ; 


i O r~r h~\ W 

* i b 


(17) 


This expression has been found to give reasonably good 
agreement with experimental data for a variety of config- 
urations. If is suggested that values of AC lfi he calculated 
from this equation for both the design under consider- 
ation and for the design for which force-test data are avail- 
able. The difference between these, values can then he 
added (w ith the proper regard for sign) to the force-test data. 

Vertical tail. The value of C t determined from foree- 

v tni i 

test data on a similar design can he corrected as follows to 
obtain (\ for the design under consideration: 


(<\ ) 

\ p t ail / design 


ro- 


( rv *- Cl 

C h L 


(18) 


Th<* results of reference .‘$5 indicate that (\ can also be 

h « i l 

affected by the location of the horizontal tail with respect to 
the vertical tail. II the two designs have approximately the 
same horizontal tail size and location, how r ever, this effect 
can he neglected . 

1 he value of ( V for a V- 1 ail can he estimated from the 

p t nil 

following empirical formula: 



where b v . tttil is the developed (not projected) span of the 
V-t ail, z v . tail is the vertical distance from the center of 
gravity to the chord of the V-tail (positive up), and V is 
the dihedral angle of tin* V-tail. More information on 

V-t ails can he found in references 30, 62, and 63. 

In the cast* of a vertical fail located on tin* wing, there is, 
in addition to tin* incremental produced by the tail 

lateral force, an incremental produced by the interference 
effect of the vertical tail on the wing. Since this inter- 
ference effect- varies greatly with span wise and vertical 
position of the tail, it should be determined from force 
tests. Usually the interference is such that a vertical tail 
above the wing gives a negative increment of (positive 
eflect-i vo dihedral) and one below r the wing gives a positive 
increment of C lft . In general, the largest interference effects 
are obtained with vertical tails at or near the wing tips. 
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THE ROLLING DERIVATIVES Cn p , Ci p , Cr p 

f '"p 

flip wing and vertical tail are the only airplane components 
that contribute appreciably to the yawing moment due to 
lolling derivative C r „^. 1 he contributions of (lie fuselage and 

horizontal tail can usually be neglected. 

Wing. — The contribution of the wing to (7„ , can be esti- 
mated from the formula and charts of figure 10 which were 
taken from reference 81). Although these charts apply 
strictly only to wings having a taper ratio of 1 .0, experimental 
data have indicated that they will also provide fairly good 
estimates for taper ratios of 0.50, 0.25, and 0. In the esti- 
mation formula 


<7„ 


C , ( AC V) 2 




( 20 ) 


tlie value of (^7o 0 ) a should bo dotonninod, if possible, from 
force-test data obtained at high Reynolds number on the 




T-p 



b/z 

tiGi. hk 9. KlToct of dihedral angle on C\ ^ for the case of subsonic inconiprcssilile flow (Cj 
in radians; V in <iegrees). Taken from reftTences 5U and fls. ^ 


where 




(%•). 


full tpan I' 


C l 


M +4) cos A 

COS A {jd r ) Am . 0 


wing under consideration, since low Reynolds number data 
might indicate values of (C D(J ) a that are* too large. For the 
case of smooth wings with a large leading-edge radius and 
low or moderate sweep, it is suggested that (C t>0 ) a for the air- 
plane be assumed to be zero at all lift, coefficients up to the 
stall. I bis assumption will result in larger negative values 
of (J„ p than would be estimated from low Reynolds number 
data on (C„ 0 ) a and consequently should lead to conservative 
dynamic stability results since an increase in C x in the nega- 
tive direction has been found to cause a reduction in dynamic 
stability. I lie value of (fh 0 ) a for highly swept wings is often 
' cry large at high lift coefficients, especially for wings with 
rough surfaces, sharp leading edges, or triangular plan form. 
For these cases, values of ( C Do ) a determined even from low 
Reynolds number data might lead to reasonably good esti- 
mates of In all these cases, however, high-scale drag 

data should be, used whenever it is available. 

Effect of high-lift devices. -The principal effect of leading- 
edge high-lift devices is to extend to a higher lift coefficient 
the linear variation of (Jn p with lift coefficient. The formula 
and charts of figure 1 0 are directly applicable to this case. The 
effect of t railing-edge high-lift devices is not so straight- 
forward, but experimental data have indicated that the 
formula and charts of figure 10 also give reasonably good 
estimates in this case. 



t*usu of subsonic inenmiuessiblc flow. A = l.o. Token from inference Hit. 


,( AC o 


C L + 


( C "o) u 


( C o 0 ) u 
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Vertical tail. -The contribution of an isolated vertical 
tail surface to (\ can be estimated by the following approxi- 
mate formula which lias also been commonly used to esti- 
mate C n of a complete airplane: 

Ptail 


(\ = 

than 


o 1 1 C 
b b * Ptail 


(21) 


The values of C Ya should be determined from force-test 

Plait 

data as previously discussed. Instead of the geometric tail 
length lib, it will usually be better to use the effective tail 
lentrth — C„ lC Y * as determined bv force-test data. 

* Plait! Ptail 

Formula (21) then becomes 


C 

1 n „ 
Ptail 


2 



(21a) 


In the case of the conventionally located vertical tail surface, 
however, the rolling wing produces a sidewash at the tail 
which greatly alters the tail contribution to C„ p . This side- 
wash causes the values of C„ to be much more negative 

Ptail 

than is indicated by formula (21). This effect is discussed 
more fully in reference 36 in which is also presented a method 
for estimating the sidewash. Some preliminary theoretical 
studies have indicated that the effect- of the sidewash on 
C n varies considerably with tail size and tail location 

Ptail 

and to some extent with wing plan form. A comprehensive 
experimental vcrb' calion of this theory is planned but as 
vet only a few scattered checks have been obtained. For 
the ease of the conventionally located vertical tail surface, 
the following formula has been found to give estimates of 
C that art 1 in fairlv good agreement with experimental 

PtaU 


data: 


or 



This formula is based on the assumption that C„ n is zero 

"tail 

at 0° angle of attack and varies with angle of attack in the 
same manner as indicated by formula (21). Formula (22) 
or the method of reference 36 can be used satisfactorily for 
first approximations of C n for most configurations with 

conventionally located vertical tails. For more accurate 
estimates, especially for configurations having an unusual 
tail size or tail location, experimental data should he used. 

For wings of triangular plan form with vertical tails either 
directly above or above and slightly behind the wing, exper- 
imental data have indicated that neither formula (21) nor 
formula (22) gives an accurate estimate of C n hut that 

an average of the values obtained by the two formulas pro- 
vides a fairly good estimate. 

It. is obvious that these methods of estimating C„ p are 


"»;{ 4 


only approximate' and are open to question in many cases. 
Experimental and theoretical studies arc currently being 
made to provide better methods of estimating . These 

studies indicate that the sidewash from the fuselage as well 
as that from the wing should he taken into account in esti- 
mating Cn p . • When these methods become available 1 , the 

approximate methods presented herein should lx 1 discarded. 
At the present time, however, formula (22) and reference 
36 will usually provide much more accurate estimates of 
C n than formula (21) which has been in common use up 

Ptail 

until this time. 

c i 

*p 

Wing-fuselage. — Most of the rolling moment due to rolling 
(damping-in-roll derivative) <7, of an airplane is produced 
by the wing. The. effect of the fuselage can be neglected 
unless the ratio of the diameter of the fuselage to the wing 
span is relatively large (greater than about 0.3). For large 
values of this ratio, the value of C] will he smaller than that 

p 

for the wing alone by an amount that can he estimated from 
a consideration of the area and lateral center of pressure of 
the wing area included within the fuselage. (See references 
106, 111, and 115.) 

Wing. — The damping in roll of wings has been the subject 
of many experimental and theoretical investigations. (See 
references on (\ in table 11.) As a result, some methods of 
estimating Ci have been developed which have been found 
to give reasonably good agreement, with experimental results. 
The method presented in reference 81 appears to give suffi- 
ciently accurate estimates of C) for zero lift. This method 

• V 

is extended in reference 92 to permit the estimation of (\ v 
over the normal flight range of lift coefficient. Estimation 
charts and formulas from reference 92 are presented in figure 
11 . 

High-lift devices. Experimental data have indicated that 
the damping in roll of wings at low and moderate lilt coeffi- 
cients is not greatly affected by the addition of high-lift 
devices such as t railing-edge flaps, leading-edge 1 flaps, slats, 
and slots. The principal effect of such devices is to increase 
the lift coefficient at which the sharp decrease in Ci occurs. 
The* charts and formulas of figure 11 can be used to estimate 
the Ci of wings with either full-span or partial-span high-lift 
devices with fair accuracy despite the fact that the method 
is not strictly applicable to partial-span high-lift devices. 
(See reference 92.) 

Wing-tip fuel tanks. -The use of wing-tip fuel tanks 

usually increases the damping in roll of the wing. The 
experimental data of reference 94 for unswept wings indicate 
that the magnitude of the increase varies with angle of 
attack and depends upon the wing taper ratio and on the size 
and location of the tanks. Unpublished experimental data 
indicate similar effects of wing-tip tanks on sweptback wings. 
The following approximate formula for estimating the incre- 
ment in Ci produced by wing-tip tanks at low lift coefficients 
is based on the limited amount of available experimental data 
and should not be expected to yield very close quantitative 
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estimates: 

(AC ) =(C \ /Max imu m tank diameter\ „ 

V l p) tanks \ l p) tanks off y WlUg Spall )\^ T ) 

(23) 

whore, for symmetrically mounted tip tanks, 

K t ~ (> 

for tanks mounted below the wing tip or forward on the wing 

tip, 

A T ~~~ 3 

and for pylon-mounted tip tanks, 

AV- 1 

Experimental data for both unswppt and swept wings 
indicate that (A/A p ) /anA;s usually bpconips smaller with in- 
creasing angle of attack and, in some eases, actually reverses 
sign at high angles of attack so that the tanks are decreasing 
rather than increasing the damping in roll. The data of 
reference 94 can be used to obtain an approximate estimate 
of the effect of angle of attack for unswept wings. 

Tail surfaces. The contribution to (J i of conventional 
type horizontal and vertical tail surfaces is\isually very small 
and, in most cases, negligible. When an airplane rolls, the 
wing produces a rotation of flow at the tail surfaces which 
reduces the already small damping moments of the isolated 
surfaces, except- in the case of the vertical tail at high angles 
of attack where the tail center of pressure is below the center 
of gravity. 

1 he contribution of an extremely large horizontal tail to 
might not he negligible and can be estimated by multi- 
plying the value of (\* for the particular tail plan form 
obtained from the charts and formulas of figure 11 by the 

factor 0.5 -g in which the factor 0.5 is included to 

account for the rotation of flow produced by the wing. 

The contribution of an isolated vertical tail surface to C t 
is given by the following approximate formula; 

2 (ft) ^ 

As in the case of / „ this formula can be modified to provide 
an approximate correction for the effect of the wing on the 
damping in roll of conventionally located vertical tail 
surfaces; 

( ~ 2 (b) 

An analysis of this expression indicates that tin* value of 



flow. Taken from reference 92 . 


where 


(U, 


I Cl* ( 

K t A eos2 A \ 


1+2 sins 


A 


A+ 2 cos^A \ _ 1 / __ C f 2 \ 

A +4 cos A / HV C/) ■* A ) 


{<i v k+-o ( c ' p )c A -o f . 0 -a. 


(5) 


A +4 cos A 
t +4 cos A 


Ci is negligible at low and moderate angles of attack 
where zjb is positive but that it might be fairly important at 
very high angles of attack where zjb is a large* negative value. 
As in the case of experimental data indicate that, for a 
vertical tail located either directly above or above and 
slightly behind a wing of triangular plan form, the value of 
^ l Ptail can t >( ‘ estimated with better accuracy by an average 
of formulas (24) and (25) than by formula (25) alone. For 
conventional tail arrangements, however, formula (25) gives 
better correlation with experimental data. 
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Wing. — The following formula for the derivative C Yv 
(lateral force due to rolling) from reference 89 is based on 
experimental data and is the same as that presented in 
reference 25 except for an additional correction to account 
for tip suction: 


f V p _ A + eos A 
(; L ’/l +4 cos A 


tail A + 


A 


(26) 


The data of reference 89 show that this formula applies only 
for lift coefficients below that at which the drag factor 

(' I} — L begins to increase. At higher lift coeflieients the 

experimental data indicate smaller values of ( Yp than given 
by formula (26). For these cases an approximation of the 
value of (\ can be obtained from the experimental data of 
reference 8th As in the ease of , the break in the variation 
of C y with lift coefficient, should he expected to occur at 
lower lift coefficients for wings having sharp leading edges or 
rough surfaces and for wings t ested at low Reynolds numbers. 
Vertical tail. -The discussion concerning C n and 

** tail 

C, is also applicable to C Y . The value of ( y Y 

Vtatl p tail 

for an isolated tail surface is given by the formula 





(' 


Y » 


il 


(27) 


This formula can be modified as follows to account approxi- 
mately for the effects of wing sidewash in the ease of a con- 
ventionally located vertical tail: 


C Y 


2 [& (b) a . o_ 


Cr. 


(28) 


An average of formulas (27) and (28) can be used for tails 
located either directly above or above and slight ly behind the 
wing. 

THE YAWING DERIVATIVES Cn r , C/ r , Ciy 

% 

Wing-fuselage. — In the past, the contribution of the wing- 
fuselage combination to yawing moment due to yawing 
(damping in yaw) derivative (\ r has usually been found to 
be small compared with the contribution of the vertical tail. 
The fuselage contribution to the damping in yaw depends, 
of course, on the relative size of the fuselage and wing. In 
the past, the relative size of these components has generally 
been such that the fuselage contribution could be neglected. 
(See references 85 and 86.) For some recent designs which 
have a large fuselage relative to the wing, however, the 


fuselage contribution to C„ r is important. In the case of 
fuselages having flat sides or having a flattened cross section 
with the major axis vertical the fuselage contribution may 
also be important and some fuselage contribution to C„ r 
should be assumed, especially at high angles of attack. On 
the other hand, experimental data have shown that a 
flattened cross-section fuselage with the major axis hori- 
zontal can have negative damping in yaw at moderate and 
high angles of attack. 

The contribution of the wing to (\ r can be estimated from 
the formula and charts of figure 12 which were taken from 
reference 25. Values of ( V> 0 for the wing should be estimated 
from force-test data. For values of 7/c greatly different from 
zero, the charts of reference 25 can he used. The formula and 
charts of figure 12 are not considered reliable at high angles 
of at tack, especially for swept wings. The use of experimental 
data from the references on listed in table IT is recom- 
mended in this case. 

The effect of partial-span inboard flaps on (\, r can usually 
be neglected. (See reference 85.) The effect of full-span 
t railing-edge or leading-edge high-lift d evict's can be esti- 
mated satisfactorily from the formula and charts of figure 12. 
Values of C D in this case are, of course, for the wing with 
the liigh-lift device installed. 

Vertical tail. The contribution of a conventional-type 
vertical tail to C„ T at low and moderate angles of attack can 
be estimated from the formula 



or, with the effective tail length — ( 1C V* substituted 

’ PtaH; p tail 

for the geomet ric tail length i/b, 


C. 


=9 C c jO! 


(29a) 


The alternative method of estimating (L presented in 

'tail 

reference 83 will probably provide better estimates than 
formula (29) in the higher angle-of-at tack range. The experi- 
mental values for ( y n presented in reference 85 for power- 
on or flap-down configurations are 30 to 40 percent greater 
than the values predicted by formulas (29) or (29a). These 
differences are attributed to lag of sidewash effects in the 
free-oscillation tests used in measuring (7„ r . In estimations of 
C n for stability calculations, similar lag of sidewash 

T tail . , 

effects should be assumed if the oscillatory mode is of primary 
importance but no lag of sidewash should be assumed if the 
aperiodic mode is most important. 

Methods for estimating the (L for wing-tip vertical 

T tail 

tails are presented in references 72 and 85. 
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F I (jure 12 .— Concluded. 


C *r 

T1 te wing and vortical tail are the only airplane components 
that contribute appreciably to roll ing-moment-duc-to-ya wing 
derivative C t of an airplane. The contributions of the 
fuselage and horizontal tail can usually be m^lect e<l. A 
semiempirical method for estimating C if is presented in 
reference 88. This method involves the use of experimental 
data on the parameter C lj3 to correct the theoretical values of 
Ci given in reference. 25 and to estimate the value of 

7 wing 

c, 

7 wing 

Wing. — The formula of reference 88 and the charts of 
Ci /C L from reference 25 for estimating C lr are given in 

t! ■ 7 wing 

figure 13. The 4 values of C^jC L to be used in the charts can 
be obtained from figure 8. For taper ratios less than 0.25, 
values of O lr jC L and C lfi fC L for a taper ratio of 0.25 can be 
used. The value of (\ „ used in the formula should be the 

Pexp 

same as the value of (\ q estimated from experimental 
data by the method indicated in the section on C In the 
case of Ci r y however, (unlike the case of C i0 ) conservative 
dynamic stability results will usually be obtained if the 
smaller values of the derivative (based on low-scale experi- 
mental data) are used instead of the larger (theoretical) 


values. This difference is a result of the fact that either an 
increase in the normally negative value of Ci & or a decrease 
in the normally positive value of Ct r can cause reduction in 
dynamic stability. As pointed out in reference 88 the esti- 
mation procedure shown in figure 13 appears to account 
satisfactorily for the effects of high-lift devices, wing dihedral, 
and airfoil section, at least for sweptback wings. This pro- 
cedure is directly applicable to midwing configurations but 
should not be used for high- wing or low-wing configurations 
because changes in wing position produce much greater 
changes in than in C% r . Work is now being done to 
develop an estimation method for the effect of wing position 
on Ci y but until this method is published the following pro- 
cedure is recommended for estimating Ci T of high-wing and 
low-wing configurations: Adjust the value of Ci & to cor- 
respond to that of a midwing position. Then use this value 
of C i& to estimate a value of C\ r that will apply to any wing 
position. 

Vertical tail. — The contribution of the vertical tail to 
Ci is usually estimated by the formula 



where C Ya is preferably obtained from force-test data. 




25 


METHODS FOR CALCULATING LATERAL STABILITY AND ESTIMATING STABILITY DERIVATIVES 


— C n iC Y is substituted for the geometric tail length lib, 

& tail / &tail 


C y 

r tai l 


2 r. 


"“tail 


(32a) 


The discussion of lag-of-sidowash offsets for C„ apply 

r ta il 

also to Cy 

r tail 

EFFECTS OF MACH NUMBER 

The effects of Mach number on the lateral stability deriv- 
atives have been treated theoretically in many investigations 
(see table 11) but very little experimental data have been 
obtained to verify this theoretical work. Moreover, only 
a small part of this experimental work has been covered in 
published reports (reference 114) because most of it is clas- 
sified at the present time. It appears, therefore, that 
estimates of the lateral-stability derivatives for the time 
being will have to be based largely on theoretical work. 

The effects’ of Mach number on the stability derivatives 
can be usually considered negligible for all airplane com- 
ponents except the wing and vertical tail. For the low-lift- 
eoeflicient condition in the case of many high-speed airplanes, 
the vertical tail contributes more than the wing to all the 
stability derivatives except C Xp . For this reason, in calcu- 
lations for transonic or supersonic speed conditions it is 
especially important to know the effects of Mach number on 
the vertical-tail lift-curve slope or (l Y$t - 

Wing. — The effects of compressibility on the subsonic 
stability derivatives of the wing can be estimated by the 
formulas of reference 26. The values of the supersonic 
stability derivatives for some wing plan forms can be 
estimated by the references tabulated in table II. In this 
table the derivatives are grouped according to the type 
of wing plan form and to the particular derivatives covered. 
A helpful summary and discussion of the effects of Mach 
number on the derivatives for several different wing plan 
forms are presented in reference 106. A summary of the 
theoretical lift-curve slope, damping in roll, and center-of- 
pressure characteristics of various wing plan forms is pre- 
sented in reference 110. In the cases in which the theory 
shows large or abrupt changes in a stability derivative with 
changes in Mach number (for example, fig. 10 of reference 
106) special care should be taken in estimating the deriva- 
tive in that particular Mach number range. The abrupt 
changes should be smoothed or faired out in a manner similar 
to that suggested in the following section for estimating 

CV, I#W . 

In some cases, experimental data for supersonic speeds will 
be available on the sideslip derivatives and on the damping- 
in-roll derivative Ci . In such cases the experimental data 
should he used in preference to the theory. Some experi- 
mental results have indicated that the effect of the vertical 
location of the wing on the fuselage on the derivative 
might be greatly different at supersonic speeds from that at 


subsonic speeds. Since no methods are presently available 
for estimating this effect for the supersonic case, it appears 
that, at least in the case of high-wing and low-wing designs, 
force-test data are necessary for obtaining an accurate 
estimate of Ci fi . 

Vertical tail. — The sideslip derivatives produced by the 
vertical tail at transonic and supersonic speeds can be esti- 
mated theoretically hut should be obtained from force-test 
data whenever possible. These sideslip derivatives can be 
used to estimate the tail contributions to the other deriva- 
tives as pointed out previously. In estimates of the value 
of (V for transonic and supersonic speeds, corrections 

&tuil 

must be made for the efleet of Mach number on the lift- 
curve slope of the tail, and these corrections should account 
for any differences in the end-plate effect of the horizontal 
tail on the vertical tail. 

For Mach numbers below about 0.8 or 0.9 and above 
about 1.6 or 1.8 the effort of Mach number on the lift-curve 
slope of the vertical tail can be estimated satisfactorily from 
the theoretical values of references 26, 84, and 110. Since 
experimental data indicate that theoretical values of lift- 
curve slope are usually too high for Mach numbers from 
about 0.8 or 0.9 to about 1.6 or 1.8, the empirically deter- 
mined fairings shown in figure 14 are recommended for use 
as a guide in the use of the theory to obtain approximate 
estimates in this Mach number range when force-test data 
are not available. 

Experimental data have indicated that for vertical-tail 
configurations which have a tail length (distance from the 
center of gravity to the tail center of pressure) that is rela- 
tively short in terms of tail chords, the rearward shift of the 
tail center of pressure at supersonic speeds can cause an 
appreciable increase in the tail length and consequently an 
appreciable increase in the magnitude of some of the tail 
derivatives. Theoretical center-of-pressure positions for var- 
ious plan forms at supersonic speeds are given in reference 
1 10 . 

EFFECTS OF POWER 

On the basis of existing information, the effects of power 
on the lateral stability derivatives appear to be negligible 
in the case of jet-propelled airplanes but these effects are 
often very large in the case of single-engine propeller-driven 
airplanes. Methods are available for estimating some of 
these power effects but in most cases experimental data are 
necessary for making a satisfactory estimate. The eflects 
of power can he broken down into two general classes: 

(1) The effects of the lateral force produced by the pro- 

peller itself 

(2) The effects of the propeller slipstream on the wing, 

fuselage, and vortical tail of the airplane 

Effects of propeller lateral force. — A method of estimating 
the propeller-lateral-force derivative C Yfi is presented in ref- 
erence 31 which is based on the work of references 32 and 33. 
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Figure 14. — Examples of sug£i'sti><l fairing of liu*orotit;al values of lift -curvo slope for use in 
estimating values for the vertical tail in the transonie ranee. 

The contribution of Hut propeller lateral force to the other 
stability derivatives can be estimated from this derivative 
by assuming that the propeller is effectively a vertical tail 
surface and by using the expressions for the tail contribu- 
tion to the various derivatives presented in the preceding 
sections. Some experimental data on the effect of wind- 
milling propeller on all of the derivative's an* presented in 
reference G6. 

Effects of propeller slipstream.— The effects of propeller 
slipstream on the lateral-stability derivatives are usually 
much greater than the effects of propeller Literal force in 
the case of single-engine tractor airplanes. The slipstream 
effects on the wing, the fuselage, and the vertical tail can be 
considered as three independent- effects. 

The slipstream effects on the wing can usually he neg- 
lected except for the derivatives Q x and C lfm Experimental 


data showing the decrease in effective*, dihedral ( — 
with power for single-engine airplanes are presented in ref- 
erences of), 50. 57, 70, and 82. It appears highly desirable 
to determine this effect of power experimentally because 
interference effects make accurate estimations of the effect 
very difficult. The effect of the slipstream on the value of 
Ci cannot be estimated from the data on C. as 

anno . • Pwi»o 

described in the section on (\ r . In fact, this procedure 
would probably give the wrong sign for the increment of 
CV contributed by the slipstream. An approximation of 

'icing 

this increment might be obtained by estimating the slip- 
stream velocity arid the lateral displacement of the slip- 
stream caused by yawing. Usually the power effects on 
tV and C, will be greatest for the flap-extended 

configuration. 

In the ease of the single-engine airplane the effect of the 
slipstream on the fuselage is usually to increase negatively 
the values of (\ tft and C Vff . (See references 55, 56, 57, 78, 
76, and 78,) Since no accurate methods of estimating these 
slipstream effects on and C v are available, it is necessary 
to determine them from force-test data. 

The effects of tin* slipstream on the vertical tail are often 
very important and should also be determined from experi- 
mental data, if possible. The increase in dynamic pressure 
at the tail caused by the, slipstream is treated theoretically 
in reference 119 and is illustrated by the experimental data 
of references 5] , 55, 50, 57, 73, 70, and 78. Tin* experimental 
data of reference 78 also show that the propeller slipstream 
can cause a destabilizing side-wash at the tail which will tend 
to reduce tin* stabilizing effect of tin* increased dynamic 
pressure at the tail. Since these data indicate that slip- 
stream effects on the vertical tail vary greatly with airplane 
configuration and propeller arrangement (single or dual rota- 
tion), use of* experimental data appears to la* the only satis- 
factory estimation procedure at present. 

Suggested estimation procedure for power effects. — The 
following procedure is suggested for estimating power effects. 
Obtain force-test, data for tail off and tail on. Use tail-on 
data directly for (' Y ^ ( r v&7 and C t& . Estimate rolling and 
yawing derivatives as follows: 

(1) Estimate fV from reference 31 and use this 

p propeller 

derivative and proper linear dimensions to estimate other 
propeller deri vat i ves (rolling and yawing derivati\es) in the 
same manner as tail derivatives. 

(2) Subtract tail-on data from tail-off data to get values 

of TV , (\ , and a for the. power-on condition and 

p tnil 'tail 

use* those values to estimate the tail contribution to tin* 
other deri vat ives. 

(3) For tail-off values of rolling and yawing derivatives, 
use same values as for power-off for all derivatives except. 
Ci r . Estimate ( ^ as suggested in preceding section. 

(4) Add the values obtained in steps 1, 2, and 3 to get the 
rolling and yawing derivatives for the complete airplane. 

INADEQUACIES IN PRESENT INFORMATION AND METHODS 

In the course* of summarizing the estimation methods for 
the various stability derivatives, the need for much additional 
information on all the derivatives became apparent. In 
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particular, information is needed to aid in the estimation 
of the derivatives in the transonic and supersonic speed 
ranges. Additional work also needs to he done in correlating 
and analyzing existing subsonic data and in obtaining new 
experimental data for the development of seiniempirieal 
methods of estimating the subsonic derivatives without 
resort to force-test data. Another important need is for 
full-scale experimental results at all speeds lor checking both 
low-scale data and the existing methods of estimating deriva- 
tives. Details of the need for additional work along these 
lines an* discussed in tin* following sections. Studies should 
also be made to determine 1 ho conditions for which the list' 
of steady-slate stability derivatives in conventional stability 
equations is inadequate and to determine satisfactory 
methods of treating such conditions. 

TRANSONIC AN1) SUPERSONIC SPEEDS 

Additional theoretical work is needed on the estimation of 
stability derivatives in the transonic and supersonic speed 
ranges to covin’ the range of wing plan forms for all the 
derivatives. In particular, more work is needed on plan 
forms eurrently under consideration, such as wings having 
moderate sweep buck and taper. This need is illustrated by 
table 11 which indicates that very little material is available 
on the stability derivative's for such plan forms except, 
perhaps, for the derivative It appears from the table 

that this derivative and the triangular plan form have, in 
the past, received a disproportionate share of attention, 
probabl v because ot the greater ease with which they could 
be treated theoretically. 

The greatest need for work on stability derivatives at the 
present tilin' is probably in the measurement of the deriva- 
tives at transonic and supersonic speeds. Experimental 
data on wings are urgently needed for checking the theoretical 
work and for use in the development of empirical corrections 
to the theorv wherever necessary. Such corrections are 
particularly needed for fairing out abrupt variations of the 
derivatives with Mach number and for fairing through the 
Mach number range for which theory predicts infinite 
values. Examples of such discontinuities as indicated by 
theory are shown in figures S to of reference 100. Since 
experimental data obtained at supersonic* speeds on wing- 
fuselage combinations and on complete models have revealed 
interference effects that arc different from those obtained at 
subsonic speeds, it appears highly desirable to obtain at 
least a limited amount of experimental data at transonic 
and supersonic speeds to evaluate these interference effects. 
For example, investigations should be undertaken to deter- 
mine the effect of wing-fuselage interference on the derivative 
C t and the end-plate effect of the horizontal tail on the 

lift-curve slope' of the vertical tail. 

Most of the experimental data on stability derivatives at 
transonic and supersonic speeds will of necessity be obtained 
at Reynolds numbers considerably less than full-scale values 
and under test conditions which might render the results 
open to question in some cases, tull-seale checks in flight 
of the low-scale data and ot the estimation methods therefore 
appear to be desirable. Consequently the methods of meas- 
uring stability derivatives in flight now being developed by 


the Cornell Aeronautical Laboratory, the Massachusetts 
Institute of Technology, and the XACA should lx* extended 
to transonic and supersonic speeds when the 1 methods appear 
to be developed to a satisfactory degree of reliability for the 
subsonic ease. Some preliminary considerations involved 
in the use of these flight, techniques art* discussed in references 
120 to 123. 

SUBSONIC SPEEDS 

The methods presented in this report for estimating the 
stability derivatives at subsonic speeds depend either directly 
or indirectly on the list* of force-test data. These methods 
arc probably more reliable than methods which do not 
involve the use of force-test data on the particular design 
under consideration or on a similar design. Methods which 
do not rely on such data art' desirable in some east's, however, 
because tin* necessary data will not always he available. 

In the cast* of sideslip derivatives, empirical methods can 
probably lx* developed largely from existing information. 
In some cases it will be necessary to augment the existing 
information with new results since much oi the available 
force-test data were not obtained in a manner that would 
make tht* data readily usable for developing general esti- 
mation procedures. 

In the case of rolling and yawing derivatives, considerably 
less information is available than in the ease of the sideslip 
derivatives. Most of the information now available was 
obtained in the Langley stability tunnel, principally on wing 
configurations and to a limited extent on complete airplane 
models and airplane components other than the wing. 
Considerably more work is required, especially for com- 
ponents in combination, before satisfactory methods can be 
developed for estimating rolling and yawing derivatives 
without the use of force-test data on the particular design 
under consideration or on a similar design. 

In discussing the work necessary for developing new pro- 
cedures for estimating the stability derivatives without the 
use of force-test, data on tin* design under consideration or 
on a similar design, it, is useful to break the problem down 
into two parts: (1) effect- of individual components and (2) 
the effect of interference of the components on each other. 

The principal components to be considered are the fuse- 
lage, wing, vertical tail, and propeller. For the isolated 
fuselage, the main problem is the development of methods 
for the estimation of C n and then, perhaps, of C nr and 
C Y For the isolated wing, the main problem is to estimate 
the derivatives at lift coefficients above that at which 
separation begins. Such estimations can be made with 
reasonable accuracy for some of the derivatives by existing 
methods which make use of force- test data, but the develop- 
ment of methods which do not involve the use of force-test 
data will probably he very difficult . For the isolated vertical 
tail, the problem is to establish the effective tail area and 
aspect ratio from the geometry of the tail so that the lift- 
curve slope (or (7y ) of the tail can be calculated. Solutions 
to this seemingly simple problem have in the past become 
involved with interference effects so that, as yet, no reliable 
methods have been published for estimating C Yf} of the 
vertical tail from its geometry. For the isolated propellers, 
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the work that is needed at present is a systematic- check of 
existing methods of estimating the lateral force 1 on the 
propeller to determine the accuracy of these 1 methods. 

The principal interference ('fleets to lx 1 considered are 
mutual interference of the wing and fuselage, wing-fuselage 
interference on the vertical tail, horizontal-tail interference 
on tin 1 vertical tail, and propeller-slipstream interference on 
the wing, fuselage, and vertical tail. 'Hie mutual-interference 
effects of the wing and fuselage are probably important 
only for the derivatives and f,. A large amount 

of experimental data is available, for the sideslip deriv- 
atives but no procedures for estimating the interference 
effects on those derivative's have been reported. Wing- 
fuselage interference has very important effects on ( f Y(j of 
the vertical tail and consequently on all of the stability 
derivatives for some flight conditions. These effects result 
from the si d ('wash and change 1 in dynamic pressure at the 
tail which may result from sideslipping, rolling, or yawing. 
Although considerable data which show these 1 interference 
effects arc available, particularly for the 1 ease of sideslipping, 
no reliable methods exist for estimating the interference 
effects. Horizontal-tail interference also lias an important 
effect on fV nf the vertical tail for some 1 horizontal-tail 

p 

positions. Some 1 work on a limited number of configurations 


has been done* toward developing methods of estimating 
this effect but data are 1 required on more configurations 
before 1 the* generally applicable methods can be evolved. 
The* propeller slipstream can cause important effects on 
(\ & and (\ r of the wing, on and (\^ of the 1 fuselage, and 
on ( <>f till* tail (and consequently on the 1 tail contribution 
to all the derivatives). Some data are available for the 
effect of tin* slipstream on the sideslip derivatives but, 
because of the complexity of this problem, considerable 
additional data may be required before a satisfactory method 
of estimating the 1 slipstream effects can he developed. 

As mentioned in the preceding section, full-scale checks 
of low-scale data and of the 1 estimation methods are desirable. 
For the subsonic erase some of the checks can he obtained 
from large-scale wind-tunned tests but some 1 checks in full- 
se*al(‘ flight tests should also be obtained when the* various 
methods of measuring stability derivatives in flight have* 
been developed to a satisfactory degree of accuracy. 


I jA ngl e y A i: m ) na i ti c a l La no k at< > it v , 

National Advisory Committee for Aeronautics, 
Langley Field, Ya., December 13, 1350 . 


APPENDIX A 

EQUATIONS OF LATERAL MOTION 
The* dimensional equations for the 1 lateral motions of an airplane 1 are 1 


m j c 2 ( P ( P d<t> , ^ ; rfV dL d\p dL 


ix Jj2 ~ n + mk xz dt , dr dt 


dt 2 dp dt 


r — L c -~ 0 


, d q> dA d<t> j 2 d 2 \p dA d\p dA , . 
dr op dt dt 2 dr dt dv 


(Al) 

(A 2) 


dr d$ 

dp dt 


, 7 d\p dr dyp 
dt dr dt 


",7-(Lin)* + i» V y>*+m v-Y.-O 


If equations (Al) and (A2) arc divided by p V'-Sb and equation (Ad) is divided by pV 2 S, tbe equations of motion may In- 

expressed in the conventional nondimensional form in which they have generally been presented in XACA papers (for 
example, see reference 2): 


2pAV 


2 d 2 <j> 1 ,, d<t> 

ill*— 2 


dd>,cy K dV 1 ,, dt 
'• d.S + ^ KxZ d7~2 


(A4) 


j r d 2 <j) 1 , , d<j) jr d 2 \p 1 y y d\p . - 

~ mKxz d « 2 ~ 2 % lh + 2llhz ‘ d, ~2 ( "r T» n ^~ ( 

1 ,, d<f> d\J/ 1 <i\L ds 

-t ' i*-' ‘♦ +2 ' Z.~t ' Jj— ?>*+*.- . 

In order to convert these equations into u form which will reduce the number of arithmetical and algebraic steps in 



29 


METHODS FOR CALCULATING LATERAL STABILITY ANI) ESTIMATING STABILITY DERIVATIVES 

performin'; stability caleulat ions, equations (A4) are multi- 
plied by in pSb ami are written in the following form: 

(/> 2 -/„/!)«/> 0 
( A \> / )- II,, /))<!>- T (/>■— ll r I))lfl— lljfH— Hr 0 
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APPENDIX B 


APPLICATION OF THE LAPLACE TRANSFORMATION TO CALCULATING MOTIONS 


The applient ion ot tht‘ Laplace transformation to the 
calculation of the lateral motions of airplanes is presented in 
order to illustrate the development of the equations of motion 
in the form in which they an 1 presented in the present, report. 
This work is similar to that presented in references 5 and (>. 
In fact, it follows the presentation in reference 5 very closely. 
Reference 6 presents a brief explanation of the Laplace 
transformation and its application to solution of the equations 
of motion of an airplane. This report also makes reference 
to detailed explanations of the Laplace transformation. In 
cases wdiere modification of the equations presented in the 
present report are necessary, reference should be made to 
these texts for an understanding of the mathematics involved. 
Applying the Laplace transforms 

L(1)=V L(I > 4 >)=\ 4 >\— <t>(, 

A 


/■ a = (/£jX 2 — // . /; X)4>o+(X 2 — 7l,X)^o + I^2X(/>0)n-f MDlMo +R-C 
r :j = — ?/ / ,A^>od-(X — // r X)^o"h X/3 ( i y c 
Solving equat ions (B 1 ) by determinant s gives 



-/*x 

/T 

K,\ 3 -l r \ 2 


— «pX 

/“o 

X 3 — /i,X 2 


X 2 — V/;i x 

/*3 

X 2 — i/ r X 2 — — (tan 7 )X 


i-/ d X 

X 3 -/„X 2 

K,\ 3 -IX 


— -/^X 

K^-n P \ 2 

X 3 — /i r X 2 


X" — x 

— ?/;.X 2 ~ X 

X 2 — ?/rX 2 ^ (^ An 7)X 


which may be expressed as 


L((f>)^ 4>\ L(I) 2 <p)= — X<£ o (^<£)o 


a 0 X S + a i X 4 + (i 2 X 3 + a 3 X 2 + a 4 X + a 5 

* x = x 2 ( /lx 4 t B\ 3 +C\H I>\+ K) 


(B2) 


and multiplying each of the eq nations by X transforms 
equations (A5) from appendix A to 

(K 2 X 3 -« P X 2 )^+(X 3 -« r X 2 )^-» ? X^=r 2 

/ (’ \ r (\ ~\ y 

(— J/pX 2 — ~ X j<£x+|x 2 — y r X 2 — ( tan t)XJ >Ax+ 


where 

r l = (X a — f p X)<#> 0 -t-(/f[X 2 — / r X)^o + X(C<<»)o+ AiX(C^) 0 +f c 


Similarly, the expressions for <p x and /3x are 

b„\ :, + 6>X 4 + b 2 \ 3 + 6 S X 2 + J 4 X+ C 
X 2 UX 4 -| />'X 3 -r-rx 2 l-/A l-/t) 

foX J + CiX 3 -(-C 2 X J -(-<'3X+^4 

^ x= XMX*+/fX*+CX*+iDX + ff) 


(S3) 

(B4) 


where the expressions for the coefficients in equations (B 2 ) 
to (B 4 ) are given in terms of the mass and aerodynamic 
stability derivatives by equations (1) to (4) in the main body 
of this report. 

In order to obtain the actual variables from the trans- 
formed variables, an inverse Laplace transformation must he 
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applied. The expressions for </> x , and aiv llu* form 
">>x where ?/ x and r x are polynomials, the deg ree ol y x being 
higher than that of ?/ x . The inverse* transform of a fimetion 
of this type is 


L l 


/"A ^ '"(*«) -v 

VV ,fr{ a'(X,) 


(Bo) 


Tn this equation all of the roots X of v\ II an* assumed to 
he distinct. This assumption is valid for p x ; hut for <t>\ and 
fa, has two zero roots. (See equations (B2), (BT), 

an<l (B4).) The terms in the equations for <f> and \p resulting 
from the two zero roots arc* 

f/O 

^(0) i m)a (110) 

where 

ti=-~ x 2 

r\ 

From equations (Bo) and (Bh) then, the inverse* transforms 


of <£ x , ^ x , and (equations (B2), (Bd), and (B4)) are* 

0- | + A,a \ A, (B7) 

^ K\<- " h ' - : 1 >:< 0X2 i li 4 < ° x i | Il t a - (Ji8) 

^=rV' ffX '-l -f.V aXje T (V ffA; *T ■■(:&*** | f r - (B<)> 

I he* equations tor the* rolling velocity p and the vawing 
velocity r can he obtained from equations (117) and (Bx) by 
different iat ion 


* "1 j X 1 < ° 

r 


-1,) 

(Bid) 

=— ( /i,X|< " ,x j 
r 

1- /i 3 X ; ,r aX s •* />* 4 X 4 < oA ! j 


(Hi 1) 


where the* expressions for the coefficients of equations (B7) 
to (111 1 ) are given by equations (0) to IS) in the* sec. t ion 
entitled ‘ k ( ’alculat ion of Motions.” 



APPENDIX C 


SOLUTION OF BIQUADRATIC EQUATION 


Many methods are available, of course, for solving for the 
roots of a biquadratic equation. For example, there are 
Horner’s, Ferrari’s, ’Bernoulli’s, Descartes’, and Hitchcock’s 
methods; various methods of solution by trial; and also 
various graphical methods such as that given in reference 1. 
Solution by trial in which synthetic division is used, however, 
is recommended as being the simplest method for most 
lateral stability work. ’Flu 4 characteristic equation for the 
lateral motions of an airplane 

H X 4 + /fX" + ( r \ 2 \ 1 >\ + 0 

generally has two real roots and a pair of conjugate complex 
roots. For these cases the two real roots can be factored out 
easily and the remaining quadratic solved for the conjugate 
complex roots. In the few eases for which all four of tin 4 
roots of the characteristic equation arc 4 complex, Descartes’ 
method can he 4 used to factor the 4 biquadratic equation into 
two quadratics. When there arc 4 rcu-1 roots, solution by 
Descartes’ method requires more time than factoring out 
the 4 real roots singly and consequently is not recommended 
for general use. These methods of solution arc* explained 
in the 4 following sections. 

SOLUTION BY TRIAL BY MEANS OF SYNTHETIC DIVISION 

Solution for real roots by trial by means of synthetic 
division consists of successive approximations of a root and 
cheeking by synthetic division until the root is determined 
to the 4 desired degree of accuracy. This cheek by synthetic 
division is based on the fact that if a is a root of a polynomial 
/(/) then x—a is a factor of /(/) and consequently no re- 
mainder is left when f(x) is divided by x—a. 

The method of solving the 4 stability biquadratic equation 
by trial with synthetic division is explained in three steps in 
the following sections. First, the rule for synthetic division 
and a numerical example are given. Second, the specific 
use of synthetic division for factoring a biquadratic is illus- 
trated by a simplified example for which the roots are known. 
This example shows how the 4 cubic and then the quadratic- 
factors of the biquadratic arc 4 obtained. Third, the use of 
synthetic division in extracting the roots of a representative 
characteristic stability biquadratic is illustrated with special 
reference to methods of making the first approximations of 
the real roots. 

Explanation of synthetic division. — Synthetic division is 
explained in almost all algebra text books but is presented 
herein for the 4 convenience of the reader. The rule for 
synthetic division may be given as follows: 

Assume that a polynomial in x (J(x)) is to be divided by 
x—a; write the coefficients of the polynomial in order, 
supplying 0 when a coefficient is lacking. 


Multiply a by tin 4 first- coefficient, and add (algebraically) 
l lie’s product to the 4 next coefficient. 

Multiply this sum by a, add to t lit*- next- coefficient, and 
proceed until all the coefficients are used. The last sum is 
the remainder and also tin 4 value 4 of tin 4 polynomial when a 
is substituted for tin 4 variable 

For example, divide 4 ,r 4 + 3./+ u 3r 2 — .r — 6 by x — % 

] -| :h- i — fi 

+ ls-HW-t - iso :* 

i -i (>-)■ 2 H- fi2-i iso 


Use of synthetic division in factoring out roots. The use 
of synthetic division to factor out two known rational roots 
of a biquadratic equation is illustrated by the following 
simple 4 example. These two rational roots represent the 4 two 
read roots of the 4 charae , teristie* stability equation which, of 
course 4 , are 4 not normally known but can he approximated by 
the 4 method given in the 4 next section of this report. 

One factor of the 4 biquadratic is x — 1 so there 4 is no re- 
mainder when t fie 4 biquadratic is divided by the 4 root l 


1 + 3 + 3 — 1 — fi 

+1 +4+7+6 1 
1 +4+7+6 0 

Since the remainder is 0, x — 1 is one 4 factor of the biquadratic 
equation and x 3 +4x 2 + 7tf+G is another factor. Inasmuch 
as a cubic equation must have at least one real root, a second 
real root of the 4 biquadratic equation can be factored out 
of the cubic. For example a: + 2 is a factor so divide the 
cubic by the root —2. 


1 +4 + 7 + 0 


-2-4-0 


2 


1+2+3 0 

The factors of the 4 biquadratic then are x — 1, a- + 2, and 
/ 2 + 2j:+3. The quadratic factor can be solved for its 

roots by the quadratic formula. For example 


9 


X — 


2+V4-12 


= — 1 + i\ 


9 


Example of application to characteristic equation. — Rea- 
sonably accurate first approximations to the real roots of the 
characteristic equation can be obtained! from simple 4 formulas. 
Successively closer approximations can then be obtained by 
interpolating from the remainders. The 4 following example 
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illustrates the application of this method to obtaining the 
roots of the stability biquadratic. The biquadratic 

XM 10.43X 3 + 16.32XM-08.6X-9.10-0 

is of the form 

/1X 4 + /?X 3 + ”| D \ +- R — ( ) 


Since tin* coefficient R is generally much smaller than 
coefficient I) in lateral stability work, one of the real roots 
(usually the smaller of the two) is approximately equal to 
— Ejl ) or it may be more closely approximated by the 
equation 


X 


R 


!)- 


CR 

I) 


or for the particular case 


X 


68.6 — 


— 9.10 

( 1 0.32) (— 9. 1 0) 

68.6 


— 0. 1 29 


Approximating the root by synthetic division 


1 + 10. 43 +16.32 + 68. 6 — 9.1 0 


Approximation 


+ .13+ 1.36+ 2.3 + 9.10 .1284 2 

+ .13+ 1.36+ 2.3 + 9.14 .129 1 

1 + 10. 56+17. 68 + 70.9+ .04 1 


1 + 10.56+17.68 + 70.9 + 0 2 


For this root, the second approximation was determined by 
dividing the coefficient R by the fourth sum from the quotient 


Approximating the mot by synthetic division 


1: 10.56 : 

17.68 ) 70.9 

Approx limit 

9.48 - 

10.20-70.9 

-9.485 6 

- 9.49 

10.16 — 71.4 

-9.49 5 

- 9.48 - 

1 0.2;)- 

-70.4 

-9.48 4 

— 9.45 - 

10.50 — 67.9 

-9.45 3 

— 9.55 

9.64-76.8 

-9.55 2 

— 9.65 - 

8.78-85.9 

-9.65 1 

1-1 11.91 - r 

8.90- 15.0 

1 

1-1 1.01 I 

8.04- 

- 5.9 

2 

1 -J 1.11 - 

7.18+ 3.0 

3 

H 1.08 [ 

7.43+ 0.5 

4 

1 -f 1.07 f 

7.52- 0.5 

5 

1 + 1.075 + 

7.48 

0 

6 


For this large real root there is no simple method of deter- 
mining the second approximation as then' was in the ease of 
the smaller real root. The magnitude of the estimated root 
in this case is arbitrarily increased or decreased slightly from 
the first approximation. From the remainders determined 
from the first two approximations, a fairly close third ap- 
proximation can then he made. 

Factoring tin* quadratic equation obtained by setting 

X 2 + 1 .075X + 7.48 


9.10 

70.9 

This procedure generally provides a good second approxi- 
mation for the small real root. 

The cubic equation obtained by setting 

X 3 + 10.56X 2 + 17.68X + 70.9 

equal to zero is of the form 


equal to zero by use of the quadratic formula gives the fund 
two roots of the biquadratic equation 

. _ 1.075 + V 1.16 — 2 O TO 2 

X "“ 2 

- 

— — 0.538 + 2. 68t 


crX 3 + 6X 2 +cX+f/ — 0 

In most lateral-stability work, a real root of this equation 
will be approximately equal to —b or it may lx 1 more closely 
approximated by the equation 

. _ b :i + d 
b*+c 


The roots of the biquadratic equation may be checked by 
multiplying I be four factors to determine whether their 
product equals tin 1 original biquadratic 

(X — 0 . 1 2 84 ) ( X + 9 . 4 85 } ( X + 0 . 5 3 8 - j - 2 . 68 i) ( X + 0. 538 — 2 . 68 /) — 

( X 2 + 9.457 X— 1 .220) (X 2 + 1 .07X + 7.47) = 

X 4 + 10.43X 3 -) 1 6.32X 2 +08.0X— 9. 1 0 


or for the particular ease 


(1 0.56 ) 3 + 70.9 
(10. f>(>) 2 + 17.68 


— 9.65 


SOLUTION BY DESCARTES’ METHOD 

Descartes’ method of solving a biquadratic equation is 
particularly useful for solving equations which do not have 
any real roots. This method is explained in most text 
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books on advanced algebra and theory of equations. In 
general, the. method consists of reducing the biquadratic 
equation to a cubic, (‘({nation which can be solved easily. 
One root of the cubic equation is used to form two quadratic 
equations tin 1 roots of which are used to obtain the roots of 
the biquadratic equation. 

Method. -Keduce the general biquadratic equation 
. 1 X 1 + J W + ( V + /> X ■ -t • K= 0 

to the form 

X 4 b X :i - f- cX 2 -( (/X -\ -e — 0 

by dividing by A, 

Obtain the values of q, /*, and ft from the following equa- 
t ions: 




'/-y+i 63 


s=e _M + h2( ~ 3 e 

4 + 10 256 


and form the (‘({nation 



and solve this cubic equation in A for one of its roots rVO. 
Solution by trial by means of synthetic, division is recom- 
mended. Determine the values of I and m from the equa- 
tions 


m - 


i ?.+ 9 r 2 + JL 
2 ' '4./: 


Substitute' the values of l and m and the value of / used in 
obtaining / and m in the (‘({nations 


y*+2xy + l=Q 
y 1 — 2 jry + m =- 0 

and solve these quadratic (‘({nations for their roots y from 
which the roots of the biquadratic (‘({nation may be obtained 
from the following relation: 


APPENDIX D 

SPECIAL NOTATION USED IN CALCULATING MOTIONS WHEN THE CHARACTERISTIC EQUATION HAS COMPLEX ROOTS 

When two of the roots \ and X 2 are conjugate complex, the coefficients A ] and A 2t B x and Jl 2 , C\ «iul (\ will be conjugate 
complex. If R + Ii is one of tin* roots X] and if the powers of Xj are expressed as 


then 


X 1 * = /f*+ Id 
Xj = R\-\- 1 \i 


X 1 2 =J? a +/ 2 / 


X] 3 — R% ~h I :d 

xd =/f 4 H-' 1 4 l 

Xi S = /?5-( IJ 

Substitution of the root R+Ii in the expression for A x gives 

{ajl h + a Jt 4 + a 2 /f 3 + a Mi ~r u 4 


_ ^o'MT | u 7 ii 3 , s | 0 l ^ « 1 ± \ T « '2 * 3 T «3 / 2 + <* 4_/l) 1 

* 1 — ( (> /I /£ 5 + 5 # /ih + 4 + 3 A*2 + 21?/? 1 ) + ( t> /I / 5 + 5 / 4 + 4 3 + 3 // /a + 2 Z? /j ) i 


The division of these complex numbers is indicated by the equation 

x j + y l i = x , x 2 ±y_ilh , x 2 y ■ — s . 1/2 

^2+?/2i *2 2 + 2/2 2 -F2 2 +?/2 2 

It is evident from these relations that A x is a complex number. In this case new symbols are used to represent the real and 
imaginary parts of A x as follows: 

Ai=R a ^ I A i 


A* is the conjugate of A x and will be referred to as 


A <2 — R a I A 


B}" proc-edures similar to those for the A coefficients, 

( b { JR>+ bill 4 + b 2 Ii's+ bjt 2 + h A R\ + i:,) + (4o/r»+ b x / 4 + b 2 h + bA2 J rbJ l )i 

(fiAIh+oliRA A(^lh + ^im 2 ^2KR,)+(AAh+^lU 4 + A(^h^Dh + 2KI l )i 
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which mu v ho referred to as 


Hi, \ Ini 


^=11,-1,,' 

r t'\R\-\ cji\) + (<" (1 /- + r, /, } ( >I> j r-<7. i v A Ii)i 

~H\AR,+ » filial M'Hy ADH, | 2A7f,) 1 «Et A // / , ( 4('h+'M)I,-\ 'iKIji 


which may he referred to as 

(\ -Kc + Id 

and 

a=R c -f c i 

Similar analysis shows that, if the roots and A., are also conjugate compleN <| uanl it ies (X ;< = /f' + Z'/ and X 4 //' - /'/!), than 

yla = /r jl 4-/ / ,i‘# 

and 

a 4 -r' a -i' a ; 

where 

y __ + \R\ -f- (i jl'x 4^1 + «*>)d- (^7o/ , -,d ^ 77 r^Z'.-H • CI 3 / / ad-«4/ / i)?’ 

'(tL-1 /?'.,+ o/?/r,d 4r/r r -i | 2/^o+(h/l/', | 5/*/v-4r/YI-3Z>/' 2 +2Zi77)i 

Also, 

/Wf'n+zy 

and 

/#4-/?'*-ZV 

where 

yi . ( -^oZ£ s - ’' h\R 4 "f- 6:>/^ a I 63/?^ -|- A 4^Z / i 6 r,) -f- (A n /' :> d />7' 4 | Aj/^~h 6 ;1 /7 1>J\ )! 

«>AI?' :r \ f>/fZf' 4 1 4( I{', | »/>/? / 2 + 2/i'i?' 1 j-h((L47'r,-i -r,///' 4 -i 4f7' :r f 3 />/',-{- 2 /t7 

Similarly, 

ZWf'cd /'<?/ 

and 

(h-ii'c-I'd 

whore 1 

-f "fi/fVl o>tfYt ■C 3 Z ?'2 + C 4 /Z'i) + (( , 0 /'r ) d-r|/ / 4 -) r : 7 ' :i ( 

((3 A R' :> + 5 BIl\ -\- 4f 7f '* + 3 TWP7+ 2/?A' ,) +((L*1 / -\ 5 HI \ ~f- 4( 7 '* -4 3 1)I' 2 A, - 2 El ',)?* 
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